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A RAPID PERTURBATION PROCEDURE FOR DETERMINING 
NONLINEAR FLOW SOLUTIONS; APPLICATION TO 
TRANSONIC TURBOMACHINERY FLOWS 

Stephen S. Stahara, James P. Elliott, 
and John R. Spreiter* 

Nielsen Engineering & Research, Inc. 
Mountain View, CA 


SUMMARY 


An investigation was conducted to develop perturbation 
procedures and associated computational codes for determining 
nonlinear flow solutions, with the objective of establishing a 
method for minimizing computational requirements associated with 
parametric studies of transonic flows in turbomachines. The 
theoretical analysis involved the development of a rapid method 
for calculating first-order changes in nonlinear flow solutions 
due to variations of an arbitrary geometrical or flow parameter. 

The procedure developed and evaluated, referred to as the 
direct correction method, was found to be capable of determining 
highly accurate approximations to families of strongly nonlinear 
solutions which are either continuous or discontinuous, and 
which represent variations in some arbitrary parameter. The 
method consists of defining a unit perturbation by employing two 
nonlinear solutions which differ from one another by a nominal 
change in some geometric or flow parameter, and then using that 
unit perturbation to predict a family of related nonlinear 
solutions over a range of parameter variation. Coordinate 
straining is used in determining the unit perturbation to 
account for the movement of discontinuities and maxima of high- 
gradient regions due to the perturbation. While simultaneous 
multiple-parameter perturbations can be treated by the method, 
the theoretical development and results presented in this 
initial study are for the single-parameter perturbation problem. 

Although the procedure is generally applicable, the results 
reported here have been directed toward nonlinear aerodynamic 
applications. Attention is focused in particular on transonic 
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flows which are strongly supercritical and exhibit large surface 
shock movement over the parametric range studied; and on sub- 
sonic flows which display large pressure variations in the 
stagnation and peak suction pressure regions. Flows past both 
isolated airfoils and compressor cascades involving a wide 
variety of flow and geometry parameter changes are reported. 
Comparisons with the corresponding 'exact' nonlinear solutions 
indicate a remarkable accuracy and range of validity of such a 
procedure. Computational time of the method, beyond the deter- 
mination of the base solutions, is trivial. 
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1 . 


INTRODUCTION 


Given the remarkable growth in capability of advanced 
computational methods for the determination of a spectrum of 
nonlinear phenomena in such diverse disciplines as fluid 
dynamics, structures, and nuclear physics to name just a few - 
a capability which has already made many difficult calculations 
routine and which is certain to improve in the future - it is 
apparent that a need exists for complementary methods capable of 
alleviating, at least in part, the usage limitations imposed on 
these methods by their run times. The need becomes particularly 
compelling when large numbers of related cases are required as 
in parametric or design studies. Techniques such as direct 
acceleration procedures provide an important means of reducing 
computer time by improving computational efficiency of the 
solution algorithm, but these and similar methods, which enhance 
the solution algorithm itself, represent only a partial answer. 
What is most desirable is a means to minimize the actual number 
of separate calculations required in a particular application by 
extending, over some parametric range, the usefulness of each 
individual solution determined by these computationally expensive 
procedures , 

Consequently, the basic motivation underlying this study is 
to extend the usefulness of such numerical solutions computed for 
specific turbomachinery configurations and flow conditions with 
a view toward reducing the computational requirements now neces- 
sary, The nature of the present investigation is both explor- 
atory and developmental in the sense that aspects of the 
procedure such as validity, range of application, and economy 
will be investigated, and a computational code embodying all the 
results of the study will be developed. 

Two fundamental methods for accomplishing such a pertur- 
bation procedure are available; a classical approach involving 
posing and solving linear perturbation equations; and a direct 
correction method employing two or more nonlinear base solutions. 
In this report, both of these methods are discussed; and an 
evaluation of the latter method, based on a large number of 
different applications, is made. 

A crucial aspect of such perturbation methods is their 
ability to accurately treat regions where either discontinuities 
or high gradients exist. For the results presented here coordi- 
nate straining is introduced as a means of accounting properly 
for the displacement of discontinuities due to an arbitrary 
change in some solution parameter. This is shown to result 
in highly accurate perturbation predictions in the vicinity of 
the discontinuity. That idea has also been extended to improve 
predictions in the vicinity of other high-gradient regions. 
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Although the procedures developed are generally applicable, 
the specific results reported here are for aerodynamic appli- 
cations. Since one of the primary objectives of this study was 
to provide a definitive proof-of-concept of such a perturbation 
method, a large variety of perturbation results based on 
transonic small-disturbance and full potential solutions were 
studied and are presented for nonlinear subsonic and transonic 
flows past both isolated airfoils and compressor cascades. In 
order to enable a critical evaluation of the range of validity 
and accuracy of the straining procedure, emphasis was placed on 
transonic flows which are strongly supercritical and exhibit 
large surface shock movement over the parametric range studied; 
and on subsonic flows which display large pressure variations in 
the stagnation and peak suction pressure regions. 
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2 . 


ANALYSIS 


2.1 Perturbation Concept and Methods 

The basic hypothesis underlying the present procedure is 
that a range of solutions in the vicinity of a previously deter- 
mined or base solution can be calculated to first-order accuracy 
in the incremental change of the varied parameter by determining 
a linearized unit perturbation solution Q defined according to 
the relation ^ 

Q 

Approximate solution for 
conditions differing from 
those of the base solution 
by an amount characterized 
by e 


r — : \ 

Base solution 

for some flow 
quantity Q 

The effectiveness of such a method, of course, depends 
upon the ability of the relationship defined by equation (1) 
to remain accurate over a range e of practical significance, and 
the fact that the unit perturbation Qp need be determined only 
once. The significance of the unit perturbation Q is obvious. 

It represents the local rate of change of the base^flow solution 
Q with respect to the particular quantity, say q, perturbed; 
tRat is Qp = (9Q/9p)^. 

Two generic methods exist for determining Qp, each differing 
in philosophy and having its own particular strengths and weak- 
nesses. We refer to these methods simply as the linear pertur- 
bation equation method and the direct correction method. 

The linear perturbation equation method represents the 
classical approach for performing a perturbation analysis and 
proceeds by establishing and solving a linear differential 
equation for the perturbation. Although in the present applica- 
tion, we confine out interest solely to the first-order term, the 
complete procedure represents a rational approximation scheme 
capable of continuation to any order. The method proceeds by 
expanding the dependent variables in an ascending power series 
in the incremental change e of the varied parameter, inserting 
that representation into the full governing equations and then 
assembling the result into a corresponding series of linear 
equations in ascending orders in e. Higher-order solutions 


IT 

^ Linearized perturbation 
solution for a unit 
change of e (1) 


I 

i 
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in general depend on both base flow plus lower-order solutions. 
Determination of the appropriate boundary conditions is done in 
a similar fashion. 

The power of the linear perturbation equation method is that 
it requires the calculation of only one nonlinear base solution. 
With that information, any number of individual perturbations 
can then be calculated, subject to the particular governing 
linear partial differential equations and boundary conditions 
which apply. The disadvantages are that each perturbation 
problem must be posed individually, including differential 
equations and boundary conditions. Furthermore, it may be 
necessary to simplify the governing equations and boundary 
conditions to a point where they can be solved rapidly relative 
to rerunning the base flow procedure. Moreover, the perturbation 
solutions themselves may be quite sensitive to the base flow 
solutions which usually enter into the perturbation problem 
through the differential equation and sometimes through the 
boundary conditions as well. 

The fundamental limitation of the method is the restriction 
of the range over which the perturbation procedure remains valid 
to a linear one. Since this characteristic depends upon the 
local behavior of the base flow with respect to the varied 
parameter, no general statement regarding range of validity is 
possible. Typical behavior for a given class of flows must be 
ascertained by checks with the base flow procedure. Initially 
unknown at the outset of an application with this technique, 
then, are the accuracy requirements imposed on the base solution 
by the perturbation procedure and the range of parameter vari- 
ation over which the linear assumption is valid. 

For the alternative method, the perturbation solution per 
unit change of the varied parameter, Qp, is determined simply 
by differencing two nonlinear base flow solutions removed from 
one another by some nominal change of a particular flow or 
geometrical quantity. A unit perturbation solution is then 
obtained by dividing that result by the change in the perturbed 
quantity. Related solutions are determined by multiplying the 
unit perturbation by the desired parameter change and adding 
that result to the base flow solution. This simple procedure, 
however, only works directly for continuous flows for which the 
perturbation change does not alter the solution domain. For 
those perturbations which change the flow domain, coordinate 
stretching (usually obvious) is necessary to insure proper 
definition of the unit perturbation solution. Similarly, for 
discontinuous flows, coordinate straining is necessary to account 
for movement of discontinuities due to the perturbation solution. 

The attractiveness of the correction method is that it is 
not restricted to a linear variation range but rather replaces 
the nonlinear variation between two base solutions with a linear 
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fit. This de-emphasizes the dependence and sensitivity inherent 
in the linear perturbation equation method on the local rate of 
change of the base flow solution with respect to the varied 
quantity. For many applications, particularly at transonic 
speeds, the flow is highly sensitive, and the linear range of 
parameter variation can be sufficiently small to be of no 
practical use. Furthermore, other than the approximation of a 
linear fit between two nonlinear base solutions, the direct 
correction method is not restricted by further approximations 
with respect to the governing differential equations and boundary 
conditions. Rather, it retains the full character of the 
original methods used to calculate the base flow solutions. Most 
importantly, no perturbation differential equations have to be 
posed and solved, only algebraic ones. In fact, it isn't even 
necessary to know the exact form of the perturbation equation, 
only that it can be obtained by some systematic procedure and 
that the perturbations thus defined will behave in some 'generally 
appropriate' fashion so as to permit a logical perturbation 
analysis. For situations involving perturbations of physical 
parameters, such as reported here, the governing perturbation 
equations are usually transparent, or at least readily derivable. 
Finally, in applying this method it isn't necessary to work with 
primitive variables; rather the procedure can be applied directly 
to the final quantity desired. 

The primary disadvantage of this method is that two base 
solutions are required for each parameter perturbation considered. 
Furthermore, both flows must be topologically similar, i.e., 
discontinuities or other characteristic features must be present 
in both base solutions used to establish the unit perturbation. 


2.2 Previous Applications 

Detailed studies of the linear perturbation equation method 
to sensitive transonic flows, with a view toward testing the 
method as an effective tool for reducing computational require- 
ments, have not been done. The primary reason is that such 
studies quickly become overwhelming. Each perturbation problem 
must be posed individually, subject to its own particular 
governing equations and boundary conditions; and then a separate 
computational code for the perturbation established. Generally, 
the governing equations and boundary conditions of the perturba- 
tion, even though they are linear, are more involved than those 
for the base solution. Additionally, the computational and 
convergence characteristics can pose similar or additional 
problems from those of the base flow procedure. 

In an attempt to examine some of these problems for 
transonic applications in at least a preliminary fashion, an 
application of the linear perturbation equation method to 
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transonic turbomachinery flows was made in reference 1. The 
conclusions obtained from that study were that reasonable results 
could be anticipated from the method for blade geometry changes, 
such as blade thickness and angle of attack. Less satisfactory 
results were obtained for perturbation changes in overall 
quantities, such as blade spacing and free-stream Mach number, 
a result that could be anticipated a priori since such perturba- 
tions alter the basic character of the flow more rapidly. The 
most significant conclusion of that study was the demonstration 
of the primary limitation of the linear perturbation equation 
method. That is, for sensitive flows such as occur in transonic 
situations, the basic linear variation assumption fundamental 
to the technique is sufficiently restrictive that the permissible 
range of parameter variation becomes so small as to be of limited 
practical use. Some preliminary applications of the direct 
correction method, however, displayed a significantly wider 
range of perturbation solution validity, in particular for 
strongly supercritical flows when coordinate straining was 
employed to account for shock movement. 


2.3 Coordinate Straining 

The concept of employing coordinate straining to remove 
nonuniformities from perturbation solutions of nonlinear 
problems is well established and originally suggested by 
Lighthill (ref. 2) three decades ago. The basic idea of the 
technique is that a straightforward perturbation solution may 
possess the appropriate form, but not quite at the appropriate 
location. The procedure is to strain slightly the coordinates 
by expanding them as well as the dependent variables in an 
asymptotic series. It is often unnecessary to actually solve 
for the straining. It can generally be established by inspection. 
The final uniformly valid solution is then found in implicit 
form, with the strained coordinate appearing as a parameter. 

In the original applications of the method (ref. 3) , it was 
applied in the 'classical' sense; that is, series expansions of 
the dependent and independent variables in ascending powers in 
some small parameter were inserted into the full governing 
equation and boundary conditions, and the individual terms of 
the series determined. An ingenious variation in the application 
of the method was made by Pritulo (ref- 4) who demonstrated that 
if a perturbation solution in unstrained coordinates has been 
determined and found to be nonuniform, the coordinate straining 
required to render that solution uniformily valid can be found 
by employing straining directly in the known non-uniform solution, 
and then solving algebraic rather than differential equations. 

The idea of introducing strained coordinates a posteriori has 
since been applied to a variety of different problems (see ref. 3) , 
and forms the basis of the current applications. 
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The fundamental idea underlying coordinate straining as it 
relates to the application of perturbation methods to super- 
critical transonic flows is illustrated geometrically in figure 1. 
In the upper plot on the left, two typical transonic pressure 
distributions are shown for a highly supercritical flow about 
a nonlifting symmetric profile. The distributions can be 
regarded as related nonlinear flow solutions separated by a 
nominal change in some geometric or flow parameter. The shaded 
area between the solutions represents the perturbation result 
that would be obtained by directly differencing the two solutions. 
We observe that the perturbation so obtained is small everywhere 
except in the region between the two shock waves, where it is 
fully as large as the base solutions themselves. This clearly 
invalidates the perturbation technique in that region and most 
probably somewhat ahead and behind it as well. The key idea of 
a procedure for correcting this, pointed out by Nixon (refs. 5,6), 
is first to strain the coordinates of one of the two solutions 
in such a fashion that the shock waves align, as shown in the 
upper plot on the right of figure 1, and then determine the unit 
perturbation. Equivalently, this can be considered as maintaining 
the shock wave location invariant during the perturbation process, 
and assures that the unit perturbation remains small both at and 
in the vicinity of the shock wave. Obviously, shock points are 
only one of a number of characteristic high-gradient locations 
such as stagnation points, maximum suction pressure points, etc., 
in which the accuracy of the perturbation solution can degrade 
rapidly. The plots in the lower left part of the figure 1 indi- 
cate such a situation and display typical transonic pressure 
distributions which contain multiple shocks and high-gradient 
regions. Simultaneously straining at all these locations, as 
indicated in the lower right plot, serves to minimize the unit 
perturbation over the entire domain considered, and provides the 
key to maximizing the range of validity of the perturbation 
method. 


2.4 Theoretical Formulation for 
Single-Parameter Perturbations 

In order to provide the theoretical essentials of the 
correction method, consider the formulation of the procedure at 
the level of the full potential equation, as most of the results 
presented here are based on that level. We denote the operator 
L acting on the velocity potential as that which results in 
the two-dimensional full potential equation for $, i.e. 

L[^] = 0 (2) 


If we now expand the potential in terms of zero- and higher- 
order components in order to account for the variation of 
some arbitrary geometrical or flow parameter q 
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$=$ +£$_ + ... 

° ^ (3) 

q = + iq 

and then insert this into the governing equation (2) , expand 
the result, order the equations into zero- and first-order 
components, and make the obvious choice of expansion parameter 
£ = Aq, we obtain the following governing equations for the zero 
and first-order components 


L[$ 1 = 0 
o 


L, [<1, ] + ^ = 0 

11 dq o 


(4) 


Here L]^ is a linear operator whose coefficients depend on zero- 
order quantities and 9L[<l>^]/3q represents a ' forcing ' term due to 
the perturbation* Actual forms of L]_ and the 'forcing' term are 
provided in reference 1 for a variety of flow and geometry para- 
meter perturbations of a two-dimensional turbomachine, and in 
reference 7 for profile shape perturbations of an isolated 
airfoil. An important point regarding equation (4) for the 
first-order perturbation is that the equation represents a 
unit perturbation independent of the actual value of the pertur- 
bation quantity e. 


Appropriate account of the movement of discontinuities and 
maxima of high-gradient regions due to the perturbation is now 
accomplished by the introduction of strained coordinates (s,t) 
in the form 


where 


X = s + exj^(s,t) 
y = t + ey^ (s,t) 


(s, t) 
(s,t) 


N 

I 6x^x^ (s,c) 
i=l "^i 

N 

I _ (s,t) 

i=l i 


(5) 


( 6 ) 


and e6xj^, eSy^^ represents individual displacements of the N 
strained points, and xii(s,t), y-^^ (s,t) are straining functions 
associated with each of the N strained points. Introducing the 
strained coordinate equations (5) and (6) into the expansion formu- 
lation leaves the zero-order result in equation (4) unchanged, 
but results in a change of the following form for the perturbation 
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(7) 


Li[$^] 


+ 


2 O 


oq o 


Here the operators are understood to be expressed in terms of the 
strained (s,t) coordinates, and the additional operator L 2 arises 
specifically from displacement of the strained points. In 
references 6 and 7, specific expressions for L 2 are provided 
for selected perturbations involving transonic small-disturbance 
and full potential equation formulations. The primary point, 
however, with regard to perturbation equation (7) expressed in 
strained coordinates is that it remains valid as before for a 
unit perturbation and independent of e. 


In employing the correction method, equation (7) for the 
unit perturbation is solved by taking the difference between 
two solutions obtained by the full nonlinear procedure after 
appropriately straining the coordinates. If we designate the 
two solutions for some arbitrary flow quantity Q as base Qq and 
calibration respectively, of the varied parameter, we have 

for the predicted flow at some new parameter value q (ref. 8) 


where 


Q(x,y) = Qo(s,t) + ^ [Q^ (x,y) -Q^ (s, t) ] 


X = s 

+ 

e^x^ (s,t) 

y = ^ 

+ 

CoYl (s,t) 

X = s 

+ 

^ [X-S] 

0 

II 

ft 

+ 

^ [y-t] 



0 

E — ( 

0 

^C 

I 

0 


( 8 ) 


(9) 


e = q 


q 


o 


In the following section, applications of the correction procedure 
are made to predict surface properties. Also provided are the 
particular forms of the straining functions equation (6) for 
those applications. 
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2.5 Current Applications: 


Surface Pressures 


For the current applications, we have employed coordinate 
straining with the correction method to predict surface pressure 
distributions for a wide variety of single-parameter geometrical 
flow perturbations of isolated airfoils and cascades. In that 
instance where flow properties are required along some contour, 
the solutions can be represented by 


Q(x;e) ~ Q„(s) + eQ-, (s) + 

o i 

X-- s + eXj^(s) + ... 


( 10 ) 


where x is the independent variable measuring distance along the 
contour or a convenient projection of that distance, s is the 
strained coordinate, and e a small parameter representing the 
change in some flow or geometrical variable which we wish to 
vary . 


In order to determine the first-order corrections (s) , 
we require a base and calibration solution in which the calibra- 
tion solution is determined by varying an arbitrary parameter q 
by some nominal amount from the base flow value. 

In this way, the first-order correction Q]^(s) can be 
determined as 


Q_ (x) - Q (s) 

Q, (s) = 2 ( 11 ) 

q - q 

where Q is the calibration solution corresponding to changing 
the parameter q to a new value q^/ x is the strained coordi- 
nate pertaining to the calibration solution, and qc - q© 
represents the change in the q parameter from its base flow 
value. If we now desire to keep invariant during the perturba- 
tion process a total of N points corresponding to discontinuities 
or high-gradient maxima, we can represent the solution by: 

Q(x;e) = Qq(s) + eQ^(s) (12) 

where 

Q (x) - Q (s) 

Q. (s) = ^ (13) 
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X = s + 

N 

1 £ C<Sx?)*x (s) 

. ^ T C 1 1 . 

1=1 1 

(14) 

X = s + 

N 

1 £ (6x7) *x (s) 

i=i ^ 

(15) 

£ = q 

c ^c 

■ % 

(16) 

e = q - 

% 

(17) 

£ (6x?) 

C 1 

f c o. 

= (x.-x.) 

(18) 

£ (6x?) = 

= (x- - X?) 

£ 1 1 

(19) 


Here (6x^) given in equation (18) represents the displacement 
of the ith invariant point in the calibration solution from its 
base flow location due to the selected change c in the q 
parameter given by equation (16), e(6x5) given in equation (19) 
represents the predicted displacement of the ith invariant point 
from its base flow location due to the desired change e in the q 
parameter given by equation (17), and ^ij^(s) is a unit-order 
straining function having the property that 

/ 1 k = i 

( 20 ) 

0 k i 

which assures alignment of the ith invariant point between the 
base and calibration solutions. 



In addition to the single condition equation (20) on the 
straining function, it may be convenient or necessary to impose 
additional conditions at other locations along the contour. For 
example, it is usually necessary to hold invariant the end points 
along the contour, as well as to require that the straining vanish 
in a particular fashion in those locations. All of these condi- 
tions, however, do not serve to determine the straining uniquely. 
The nonuniqueness of the straining, nevertheless, can often be 
turned to advantage, either by selecting particularly simple 
classes of straining functions or by requiring the straining to 
satisfy further constraints convenient for a particular appli- 
cation. An example of the effect of employing two different 
straining functions for a strongly-supercritical flow was 
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provided in reference 6. Here we provide additional results 
demonstrating some of the limitations of various polynomial 
straining functions and provide some comparisons with piecewise- 
continuous functions. The particular classes of straining 
functions employed were continuous polynomial and linear 
piecewise-continuous. For these two classes, the functional 
forms of the straining can be compactly written. For example, 
equation (14) becomes, for continuous polynomial straining 

N-1 

X = s + I L. (s)* (xT -x°) (21) 

. ^ ^ 1 1 1 
1=2 


where L . 

1 


are Lagrangian coefficients given by 


(s) 


N 

n 

k=l 

k^i 




( 22 ) 


whereas for linear piecewise-continuous straining, x is given 
by 


N-1 

X = s + I 
i=2 


o 

o o 

X . , T - X . 
1+1 1 



+ 





f C o \ 

(^i+1 -Xi+l) 


H (x?_|^l - s) "H (s - X?) (23) 


where H denotes the Heaviside step function. As discussed 
above, it is usually necessary to hold invariant both of the 
end points along the contour in addition to the points corre- 
sponding to discontinuities or high-gradient maxima. Consequently, 
for the results reported here, the array of invariant points in 
the base and calibration solutions have been taken as 


0 

X . 

1 


c 

X . 
1 


{0, 

o 

Xi, 

o 

^2' • ' 

o 

. . . , x^. 

1} 

{0, 

c 

Xi, 

c 

X2' •' 

c 

1} 


(24) 


where the contour length has been normalized to unity. Figure 2 
provides a summary of the various combinations of flows and 
straining functions employed. 
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3. 


RESULTS 


One of the primary objectives of the present investigation 
is to explore the accuracy and range of validity of such pertur- 
bation procedures to determine to what extent they are capable 
of providing results useful in an engineering analysis. To 
this end, we have tested the correction method with coordinate 
straining over a wide variety of different geometrical and flow 
condition perturbations, including applications to both isolated 
airfoils and compressor cascades. In particular, since the 
ability of the method to account accurately for the movement of 
discontinuities and maxima of high-gradient but continuous 
regions is essential if such procedures are to be of general use, 
emphasis was placed on transonic flows which are strongly super- 
critical and exhibit large surface shock movement over the para- 
metric range studied. Base flow theoretical solutions were 
determined from small-disturbance transonic potential (ref. 9) 
and full potential solutions (refs. 10, 11, 12). In the 
results to follow, which were selected as typical from systematic 
calculations of a much larger number of cases, the choice of 
base and calibration solutions was often made at the limits of 
validity of the procedure to observe how well the method works 
under such conditions. 


3.1 Perturbation Results for Supercritical Single- 
Shock Flows and Subcritical Flows 

3.1.1 Su percritical applications .- In figure 3, V7e present 
results for a thickness-ratio perturbation of strongly super- 
critical flows past a nonlifting cascade of biconvex profiles 
at M^ = 0.80 having a spacing-to-chord ratio of H/C = 1.0. The 
dotted and dashed results on the figure represent the base and 
calibration surface pressure distribution for t = (0.075, 0,065), 
respectively, and were obtained by solving the transonic small- 
disturbance potential equation using the code TSFOIL (ref. 9) . 

An x-grid having 48 points on the blade profile was used. These 
solutions were then used to determine the unit perturbation. The 
open circles represent the perturbation solution for t= 0.073 in 
the plot on the left and for t =0,070 in the plot on the right. 
Those perturbation results are meant to be compared with the 
solid lines in the plots which are the corresponding nonlinear 
solutions obtained by rerunning TSFOIL at the new thickness 
ratios. Quadratic straining was used with shock point and leading 
and trailing edges held invariant. The base and calibration 
flow shock-point locations for this example, as well as for all 
of the supercritical cases presented here, were determined as 
the point where the pressure coefficient passed through critical 
with compressive gradient. 
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with regard to the results, several points are noteworthy. 
Selection of a cascade rather than an isolated airfoil provides a 
more sensitive transonic flow situation. Additionally, the choice 
of a highly supercritical base and almost subcritical calibration 
solution provides both an example of extreme separation between 
the two nonlinear solutions used to define the unit perturbation, 
as well as a situation where one solution is near the limits of 
validity of the perturbation analysis. Recall that both solutions 
must be topographically similar, i.e., must contain the same number 
of discontinuities (shocks) and other characteristic features. 

We note that comparisons of the perturbation results with 
the nonlinear calculations are very satisfactory for both 
thickness ratios, with the only discrepancy being a slight disa- 
greement at the lower thickness ratio (t =0.070) at several 
points in the post-shock region. Additional calculations not 
presented here in which a more reasonable choice of calibration 
solution is made, say at t =0.070, removes that discrepancy as 
well. The main point provided by the results of figure 3 is 
that for certain classes of supercritical flows even widely 
separated base solutions can be used to provide reasonable 
perturbation predictions. 

In figure 4, we provide similar strongly supercritical 
results again for interpolation-only perturbation solutions, but 
in this instance on a somewhat finer grid. These results 
employed full potential base solutions (ref. 10) , and represent 
thickness ratio perturbations of nonlifting symmetric free-air 
flows past NACA four-digit thickness-only airfoils at M^ = 0.820. 

The body-fitted mesh employed had 75 points on both upper and 
lower surfaces, which is half again as many as in the preceding 
example. For the base and calibration flows, the thickness 
ratios were t =0.120 and 0.080, respectively. Comparisons between 
the perturbation predictions and the full nonlinear calculation 
are exhibited in figure 4 for x= 0.110, 0.105, 0.100, and 0.095. 

We note that the comparisons are remarkably good, in particular, 
in the region of the shock. The first-order perturbation 
accurately predicts both shock location and the post-shock 
expansion behavior. Reference to the coarser grid results given 
in figure 3 indicates that the finer grid resolution clearly 
enhances the perturbation result, indicating that better 
accuracy and a larger range of validity of the perturbation 
solutions can be anticipated when fine-grid base solutions are 
used to define the unit perturbation. 

In the two preceding examples, perturbation results were 
provided for interpolation-only between widely spaced base and 
calibration solutions. In figure 5, we provide similar strongly 
supercritical thickness-ratio perturbation results for extreme 
solution extrapolation using very closely spaced base and 
calibration solutions (ref. 10). The upper plots display results 
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for extrapolation downward from base and calibration flows past 
nonlifting NACA OOXX profiles with t = 0.115 and 0.120 at 
M^= 0.820. Perturbation predictions are shown for t= 0.105 and 
0.100, which represent At excursions from the base flow 
(t = 0.115) that are two and three times the parameter change 
between the base and calibration solutions (At =0.005) used to 
define the unit perturbation. For these results, comparisons 
with the full nonlinear calculations are very good. The lower 
plots display similar results for extreme extrapolation upward 
from base and calibration solutions have t = 0.095 and 0.090. 
Perturbation predictions are shown for t= 0.105 and 0.110, which 
again represent excursions from the base flow that are two and 
three times the parameter change between the base and calibra- 
tion solutions. In this instance, while comparisons of the 
perturbation results and the full nonlinear solutions for both 
cases are good, the results at t = 0.110 are beginning to display 
some not surprising discrepancies near the shock wave, indicating 
that the perturbation result is nearing the limit of its range of 
validity for this particular choice of base and calibration flows. 

The results indicated in figure 5, however, clearly demon- 
strate that not only is accurate solution extrapolation possible, 
but that for some situations even closely spaced nonlinear solu- 
tions can be used to cover a wide range of related solutions. 
Additionally, the range of parameter variation in this example 
over which the perturbation results remain accurate - i.e., 
parameter changes three times the difference between the two 
nonlinear solutions used to define the unit perturbation - is 
remarkable, and far beyond what one would anticipate for a 
first-order correction. 

Perturbation results using a more reasonable choice of base 
and calibration solutions are provided in figure 6. Those 
results involve Mach number perturbations of highly supercri- 
tical full potential (ref. 10) flows past a NACA 0012 airfoil at 

a= 0°. The base and calibration results are for M = 0.800 and 

00 

0.820, and the comparisons indicated are for perturbation results 
interpolated to M^ = 0.810 and extrapolated downward to M^ = 0.790. 
As in the case of the geometric perturbations given in figures 4 
and 5, these perturbation results are also in very good agreement 
with the nonlinear calculations at the new Mach numbers. For 
this perturbation, as well as for a number of other Mach number 
perturbations, we have separately determined the perturbation 
result in two different ways. First, we have taken cognizance 
of the fact that a Mach number perturbation alters the governing 
differential equation for the first-order perturbation from that 
of other geometric or flow parameter changes; and have used the 
suggestion of reference 6 to consider such perturbations via a 
transonic small-disturbance approximation, whereby the same 
perturbation equation can be preserved by employing a modified 
expansion parameter s . TVn alternative procedure is to treat a 
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Mach perturbation directly and interpret e as the difference in 
Mach number. We have done these calculations and compared the 
perturbation results for a number of cases using both full 
potential solutions, as for the results shown in figure 6, and 
transonic small-disturbance solutions, and have observed no 
essential difference between the two sets of results. The 
perturbation results presented in figure 6 correspond to those 
for e equal to the difference in Mach number. 

All of the supercritical perturbation results presented in 
figures 3 to 6 have been for symmetric flows and have employed 
a quadratic straining function. In figure 7, we present results 
for an angle of attack perturbation of lifting flows past a 
NACA 0012 profile at M^ = 0.70. The full potential (ref. 10) 
base and calibration solutions are at a =3.0° and 4.0°, with 
comparisons of the perturbation and full nonlinear results shown 
for a =3.5° and 2.5°. Cubic straining has been used with the 
invariant points corresponding to the lower trailing edge, stag- 
nation point, shock point, and the upper trailing edge (see 
fig. 2). We note that a =3.5°, the perturbation results are very 
good everywhere, in particular, in the vicinity of the shock and 
stagnation regions. At a =2.5°, the perturbation results are still 
very good in the shock and stagnation regions and on most of the 
upper and lower surface, but near the trailing edge a discrepancy 
has occurred. The cause of this discrepancy lies solely with the 
cubic straining function used. It is due to the fact that 
although the straining vanishes identically at the trailing edge, 
for the particular choice of base and calibration solutions in 
this example, the straining in the near vicinity of the trailing 
edge becomes sufficiently large to introduce a misalignment in the 
unit perturbation in that high-gradient region. The correction to 
this is discussed in the section describing piecewise-continuous 
straining functions. 


3.1.2 Subcritical applications. - Although supercritical 
flows are clearly of central concern in any transonic analysis 
for which the perturbation methods presented here would be used, 
applications to subcritical nonlinear flows are also of signifi- 
cance. To this end, we have applied these same techniques to a 
variety of subcritical flows to examine their accuracy and range 
of validity for such applications. 

In figure 8, we present some summary results for four 
different subcritical perturbation applications to an isolated 
airfoil. All of these results are based on full potential 
solutions (ref. 10) with quadratic straining holding invariant 
the stagnation point and the trailing edge points. The plot on 
the upper left displays comparisons for a camber line pertur- 
bation of a lifting flow with = 0.50 and a = 2° past an airfoil 
having a NACA 0012 thickness distribution and a parabolic-arc 
camber line having the maximum camber located at midchord. Base 
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and calibration flows with camber ratio h/c = 0.02 and 0.01 were 
used to extrapolate perturbation results to h/c = 0,05. Compar- 
isons with the full result is essentially exact. The plot on the 
upper right provides similar results for a thickness-ratio pertur- 
bation of a lifting flow with M^ = 0.50 and a = 2 . 0® past NACA 
OOXX thickness-only airfoils. Base and calibration flows with 
T =0.12 and 0.04 were used to provide interpolation results at 
T =0.08. Again, the agreement is essentially exact even in the 
peak suction pressure region. The plot on the lower left provides 
angle-of-attack perturbation results for M^ = 0.50 flow past a 
NACA 0012 airfoil, using base/calibration results for a =4.0°, 

2.0° to predict results at a =3.0°, with the agreement again 
being quite good. The final comparisons given in the plot on 
the lower left are for a Mach number perturbation of a lifting 
flow at a = 2° past an airfoil having a NACA 0012 thickness 
distribution and a parabolic-arc camber line with camber ratio 
h/c=0.03 at midchord. Base/calibration results at M^=0.40, 

0.60 were used to predict results at M^ = 0.55, with good agree- 
ment with the full nonlinear calculation. 

In figure 9, we present similar summary results for subcrit- 
ical perturbation applications to a compressor cascade having 
a 4% biconvex thickness distribution and a 1% parabolic-arc 
camber line blade, a pitch of t/c=0-37, and oncoming Mach 
number M^= 0.770. These results are based on the full potential 
solution procedure of reference 11 and have also used quadratic 
straining to hold the trailing edge points and stagnation point 
invariant. The plots in the upper part of the figure represent 
an inflow angle perturbation, with base/calibration inflow 
angles =47.8°, 49.8° used to predict extrapolation results 
in the plot on the left for 3i=48.8° and interpolation results 
in the plot on the right for 3j_ = 48.8°. In the lower left 
plot, interpolation results are displayed for an outflow 
angle perturbation with base/calibration outflow angles 
3 q = 31.5°, 39.5° used to predict the flow at 3o = ^5.5°, The 
lower right plot provides interpolation results for a rotational 
speed perturbation with base/calibration rotational speeds 
CO = 967,667 rad/sec used to predict the flow at co = 82 7 rad/sec. 

In all of these results, the perturbation results are good, 
including the regions near the leading and trailing edge where 
a peaky behavior due to local grid resolution is observed. 


3.2 Comparison of Continuous and Piecewise-Continuous 
Straining Function Perturbation Results 

The results presented in figures 10 to 13 illustrate the 
effect of using different straining functions to determine the 
perturbation results. Comparisons are provided for several 
strongly supercritical flows, demonstrating the differences in 
perturbation solutions between using quadratic and cubic straining 
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functions and corresponding piecewise-continuous straining 
functions . 

Figure 10 displays a comparison for a symmetric supercrit- 
ical thickness-ratio perturbation at t =0,110 for which results 
based on quadratic straining were given in figure 4, In that 
figure the open circles denote the previously obtained perturba- 
tion results using quadratic straining, while the asterisks denote 
the corresponding result when using linear piecewise-continuous 
straining. The points held invariant are the leading and trailing 
edges and the shock point. For this case there is virtually exact 
agreement everywhere between the two perturbation results as well 
as with the nonlinear result. An analogous comparison with a 
cubic straining result is provided in figure 11 where the invariant 
points are the lower trailing edge, stagnation point, shock point, 
and upper trailing edge. Displayed in that figure as open circles 
are the cubic-straining supercritical angle-of-attack perturbation 
results at a =2.5° which were previously given in figure 7. Aster- 
isks denote the corresponding linear piecewise-continuous strain- 
ing perturbation result. We note that the discrepancy near the 
trailing edge caused by the cubic straining has been effectively 
removed in the piecewise-continuous result. Moreover, the good 
agreement with the full nonlinear result which the cubic result 
displayed near the shock and stagnation regions, as well as over 
the remainder of the airfoil surface, is also obtained with the 
piecewise-continuous result. 

Finally, we have found that when employing quadratic, cubic, 
and higher-order polynomials as straining functions, for certain 
combinations of base flow shock location and shock movement 
between base and calibration solutions, particularly when large 
shock movements are involved, the polynomial straining functions 
will strain some points off the airfoil surface. This, of course, 
invalidates the determination of the unit perturbation, and 
requires that a different straining function be employed. Piece- 
wise-continuous straining functions provide a simple means of 
avoiding such difficulties. 

In figures 12 and 13, we have provided examples illustrating 
this effect for both quadratic and cubic straining functions. 

Figure 12 provides a comparison of perturbation results obtained 
using quadratic (open circles) and linear piecewise-continuous 
(asterisks) straining applied to a supercritical Mach number 
perturbation for symmetric nonlifting flow past a NACA 0012 
airfoil. Widely separated base/calibration flows (ref. 10) at 
M^^ = 0.820 and 0.750 were used to predict the flow at M^ = 0.810. 

The spurious behavior near the leading edge displayed by the open 
circles is due to the quadratic function moving points in the 
strained calibration solution off the airfoil surface. The 
piecewise-continuous results indicated by the asterisks display 
a smooth variation in that region, and provide good agreement 
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everywhere with the full nonlinear result. Figure 13 provides 
a corresponding comparison for cubic straining. Angle-of-attack 
perturbation results at M^ = 0.70 for flow past a NACA 0012 
profile using base/calibration results (ref. 10) at a =2.25° and 
4.00° are used to predict the flow at a =3.25°. The unusual 
results displayed by the open symbols near the trailing edge 
indicate that the cubic function has strained points off the 
airfoil surface in that region. However, the linear piecewise- 
continuous result corrects that problem and displays good agree- 
ment with the nonlinear calculation in that region as well as at 
the shock and stagnation point. 


3,3 Perturbation Applications to Complex 
Supercritical Flows 

In order to provide a seveire test of the perturbation 
procedure, we have applied the method to a number of transonic 
flows that are characterized by surface pressure distributions 
having multiple shock and/or high-gradient locations, such as 
those typified schematically in the lower plots of figure 1. 
Demonstration of the ability of the perturbation method to 
predict accurately such classes of flows, which are typical of 
those encountered in certain transonic turbomachinery applications, 
is crucial to the present study. In order to accomplish such 
a demonstration, we have investigated two separate classes of 
sensitive supercritical transonic flows, i.e. those with multiple- 
shock waves, and those having a single shock together with 
multiple high-gradient regions. Examples of perturbation 
results for such flows are provided below. 


3.3.1 Multi-Shock Supercritical Flows .- In figure 14, we 
present results for an angle-of-attack perturbation of super- 
critical lifting flows past a NACA 0012 profile at M^ = 0.80. 

These highly sensitive flows exhibit two shocks, one on each the 
upper and lower surface. The full potential (ref. 10) base and 
calibration flows employed are at a =0.50° and 0.20°, with 
comparisons of the perturbation and full nonlinear results shown 
for a=0.0°, 0.1°, 0.4°, and 0.6°. Piecewise-continuous linear 
straining has been used with the invariant points corresponding to 
the lower trailing edge, lower surface shock point, stagnation 
point, upper surface shock point, and upper trailing edge (see 
fig. 2) . We note that the symmetrical extrapolation result at 
a =0.0° is separately predicted from both the upper surface and 
lower surface pressure distributions, and, as can be seen, the 
results are quite good. The remaining results at a = 0.1°, 0.4°, 
and 0.6°, which represent both extrapolation and interpolation 
from the base and calibration flows, are in excellent agreement 
with the full nonlinear result. As an indication of the sensi- 
tivity of these flows, we have found that the lower surface shock 


21 


disappears at an angle of attack of approximately 0.8®; yet the 
lower surface pressure distribution is well predicted by the 
perturbation result over the parametric range studied. 

3.3.2 Supercritical Compressor Cascade Flows .- As an 
example of the ability of the method to predict a complex super- 
critical flow, in figure 15 we provide results for oncoming 
Mach number perturbation of supercritical flows past a cascade 
composed of Jose Sanz (ref. 12) profiles. For these results, 
the oncoming and exit flow angles are 30.81° and 0.09°, respec- 
tively, the blade twist is 9.33°, while the gap to chord ratio 
is 1.028. The full potential (ref. 12) base and calibration flow 
oncoming Mach numbers are M^=0.77 and 0.81, with comparisons of 
perturbation and full nonlinear results shown at M^=0.75, 0.79, 
0.89, and 0.83. Piecewise-continuous linear straining was 
employed with invariant points at the lower trailing-edge , stag- 
nation point, shock point and upper trailing edge. As with the 
multiple-shock example shown in figure 14, we note that the 
perturbation predictions are in excellent agreement with the 
-nonlinear results. In particular, we note that the perturbation 
procedure captures the variation of the plateau-like pressure 
distribution on the upper surface near the leading edge, the 
location and strength of the shock, the post-shock expansion 
region, the rapid expansion near the trailing edge, and the 
expansion on the lower surface near the stagnation point, indi- 
cating a capability for treating very general flow situations. 
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4 . 


CONCLUSIONS AND RECOMMENDATIONS 


An evaluation has been made of a perturbation procedure for 
determining highly accurate approximations to families of non- 
linear solutions which are either continuous or discontinuous, 
and which represent variations in some arbitrary parameter. The 
procedure employs a unit perturbation, determined from two 
nonlinear solutions which differ from one another by a 
nominal change in some geometric or flow parameter, to predict 
a family or related nonlinear solutions. Coordinate straining 
is used in determining the unit perturbation in order to account 
properly for the motion of discontinuities and maxima of high- 
gradient regions. Extensive perturbation calculations based 
on full potential nonlinear solutions have been carried out. 

These calculations cover a variety of flow and goemetric parameter 
perturbations involving isolated airfoils and compressor cascades 
at both subsonic and transonic flow conditions. Particular 
emphasis was placed on supercritical transonic flows which 
exhibit large surface shock movements over the parameter range 
studied; and on subsonic flows which display large pressure 
variations in the stagnation and peak suction pressure regions. 
Perturbation results for single-parameter perturbations, 
characterized by both extreme solution interpolation using 
widely separated base flow solutions and extreme solution 
extrapolation using closely spaced based flow solutions, were 
obtained in order to determine the accuracy and range of validity 
of the method. Additionally, calculation of perturbation results 
were made to investigate the effectiveness of employing piece- 
wise-continuous straining functions rather than polynomial 
(quadratic, cubic, quartic) functions. Multi-shock and other 
complex flow situations were studied in order to examine the 
capability for treating general transonic flows. 

Comparisons of the perturbation results with the corre- 
sponding 'exact' nonlinear solutions indicate a remarkable accu- 
racy and range of validity of the perturbation method across 
the spectrum of examples reported. Geometry and flow parameter 
perturbations are treatable with equal success. Solution inter- 
polation and extrapolation are both feasible. Results evaluating 
the polynomial and piecewise-continuous straining functions 
indicate that the piecewise-continuous functions are superior. 

The latter class of straining functions eliminate both the problem 
of unwanted straining in the domain of interest, as well as the 
problem of spurious straining out of the domain. Finally, it 
was demonstrated that this procedure can successfully treat flows 
containing multiple shocks and high-gradient regions by simul- 
taneously straining all of these characteristic points. 
Computational time of the method, beyond the determination of the 
base solutions, is trivial, A code encompassing these developments 
has been written for the single-parameter perturbation problem 


23 



and is included as part of this report. Based on these results, 
we conclude that such a perturbation procedure can provide a 
means for substantially reducing computational requirements in 
design studies or other applications where large numbers of 
related nonlinear solutions are needed. Further development is 
needed, however, to provide a computational tool of wide utility. 
Because of the practical need in design or parametric studies to 
consider variations in several parameters simultaneously, we 
suggest the development of the capability for multiple-parameter 
perturbations, making full use of the current developments of the 
single parameter procedure. That procedure should incorporate 
a limiting-parameter calculation whereby the parameter bounds 
with respect to each varied parameter are determined. Finally, 
in order to demonstrate their ultimate power and utility, these 
procedures should now be tested by actual application to a 
practical problem which involves the high-frequency use of expen- 
sive computational codes in order to determine a large number of 
related flow solutions. We suggest transonic turbomachinery 
blade design optimization studies as both feasible and of high 
current importance. 
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APPENDIX A 


USER'S MANUAL FOR COMPUTER PROGRAM PERTURB 


A.l INTRODUCTION 

The purpose of this appendix is to describe the operation 
of the computer code which was developed in conjunction with the 
theoretical work presented in this report, and to provide 
sufficient detail to permit convenient use and change of the 
program. The program computes and plots an arbitrary flow 
variable on a contour surface by employing the strained-coordinate 
perturbation method previously discussed. The plot package 
included in this version refers to system routines at the Stanford 
University Center for Information Processing facility. In general, 
the plotting software must be supplied by the user according to the 
requirements of his operating system. This can be accomplished 
directly by replacing or modifying the subroutines PLOT, LIMITS, 
and ROUND. 

A description of the general operating procedure of the 
program is given, together with complete description of both 
input and output. The program is written in FORTRAN IV and has 
been developed on an IBM 3033 computer. Typical run times are 
1 to 3 seconds. The storage requirements are SOK^^q. 
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A. 2 PROGRAM DESCRIPTION 


The program calculates both continuous and discontinuous 
nonlinear perturbation solutions which represent a single- 
parameter change in either geometry or flow conditions by 
employing a strained-coordinate procedure. The method utilizes 
a unit perturbation, determined from two previously calculated 
solutions ('base' and ' calibration ' solutions) obtained from 
an 'expensive' computational procedure and displaced from one 
another by some reasonable change in geometry or flow variable, 
to predict new nonlinear solutions over a range of parameter 
variation . 

This version of the procedure is configured to predict and 
plot an arbitrary flow variable (e.g., pressure coefficient) on 
the surface of a blade or airfoil , and can account for the motion 
of ; 

1. one or more critical points (shock points) , 

2. a stagnation point, 

3. a maximum-suction-pressure point, 

or simultaneously for any combination of these. 

The program is also configured to compare the perturbation- 
predicted solutions with the corresponding 'exact' solutions 
obtained by employing the same 'expensive' computational procedure 
used to determine the base and calibration solutions. 

The coordinate straining employed is piecewise linear with 
the end points and up to six interior points held invariant. At 
the option of the user, these additional interior points may be 
arbitrarily preselected, or chosen from among the minimum, maxi- 
mum, and critical points automatically located by the program 
itself . 

Critical or shock points are located on the basis of a user- 
supplied statement function defining the critical value of the 
dependent variable as a function of some single flow variable. 

The program default is with dependent variable y defined as 
pressure coefficient, with the independent variable being Mach 
number. In this case, the critical value is defined as 
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where y is the ratio of specific heats. If instead of surface 
pressure coefficient, the surface velocity distribution were used, 
then the value of y would be given by 
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Data for base, calibration, and comparison solutions (if 
available) are input as an array x(I) of coordinates and a cor- 
responding array y(I) giving the dependent variable at each 
coordinate location, where 1 < I < N and N < 200. 


Oncoming 

flow 



The leading edge is at x = 0; the data are read in beginning on 
the lower surface at the point farthest from the leading edge and 
proceeding clockwise around the surface as shown in the sketch. 
Data for the different solutions need not correspond to identical 
locations on the surface, except for the initial and final points, 
i.e., x(l) and x(N) must be the same for all cases. The program 
normalizes the x coordinates CO ^ x < 1) such that x = 0 corres- 
ponds to I = 1 and X = 1 to I = N. 
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The base and calibration solutions are searched for minimum, 
maximum, and critical points, e.g.. 



flow 


Note that the sign of dy/dx in physical coordinates is used in 
determining the critical points. For example, both critical 
points indicated on the above figure correspond to dy/dx > 0 in 
physical coordinates, since at point #1 the physical coordinate 
increases in the direction from right to left, whereas at point #2 
it increases from left to right. 

The points to be held invariant in straining are either 
selected from among those located by the program or individually 
specified by the user, after which the unit coordinate straining 
and unit perturbation are computed. 

Data for the test cases is then read in and nonlinear per- 
turbation solutions constructed from the unit perturbation. 
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A. 4 DICTIONARY OF INPUT VARIABLES 


A 


B 


LCHEK 


LECHO 


LOCO (I) 


LOCI (I) 


LPERT 


LSELCT CD 


Scaling parameter in straining procedure. A - -xCl)/ 
where xCl) is location of first data point on lower 
surface (see PROGRAM DESCRIPTION) . 

Scaling parameter in straining procedure. B = x(N), 
where x(N) is location of last data point on upper 
surface (see PROGRAM DESCRIPTION) . 

Specifies whether or not perturbation solution is to be 
compared with an exact solution. 

LCREK = 0 ... no comparison 
LCHEK =1 ... comparison 

Controls whether or not input deck is printed. 

LECHO = 0 ... no print 
LECHO = 1 . . . print 

Array of length 6 containing subscripts of user-specified 
invariant points in base solution; operational only when 
LSPEC = 1. 

Array of length 6 containing subscripts of user-specified 
invariant points in calibration solution; operational 
only when LSPEC = 1. 

Specifies type of perturbation; operational only when 
LCHEK = 1 and only affects output from plot subroutine. 

LPERT =1 ... thickness-ratio perturbation 
LPERT =2 ... angle-of-attack perturbation 
LPERT =3 ... Mach-number perturbation 

Array of length 6 of which NSELCT elements are read in; 
operational only when LSPEC = 0, and specifies nature 
of points to be held invariant according to the code: 

1 ... minimum point held invariant 

2 ... maximum point held invariant 

3 ... 1st critical point held invariant 

4 ... 2nd critical point held invariant 

5 ... 3rd critical point held invariant 

6 ... 4th critical point held invariant 

Note that critical point ordering is determined from 
order of occurrence starting at the lower surface at the 
point furthest from the leading edge and proceeding clock- 
wise around the surface (see PROGRAM DESCRIPTION) . 
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LSPEC 

LUNIT 

MO, Ml, M2 
N 

NAME 

NCASE 

NSELCT 

Q0,Q1,Q2 

TITLE 


Note that the code nujnbers can be assigned in any 
order, e.g,, 

LSELCTCD =1 LSELCTCD = 4 

LSELCTC2) = 3 and LSELCTC2) = 1 

LSELCTC3) = 4 LSELCT(3) = 3 

are equivalent, both corresponding to NSELCT — 3, with 
the minimum, and first and second critical points held 
invariant. 

Controls how invariant points in straining are 
specified . 

LSPEC =0 ... invariant points selected from 

among those located by the program, 
using the array LSELCT(I) 

LSPEC =1 ... invariant points preselected by user, 
using the arrays LOCO (I), LOCI (I) 

Controls whether or not unit coordinate straining and 
unit perturbation are printed. 

LUNIT = 0 ... no print 
LUNIT =1 ... print 

Oncoming Mach numbers in base, calibration, and 
perturbation solutions. 

Number of locations for which data are input for base, 
calibration, and comparison solutions. 

Character string of length 2 which symbolizes dependent 
variable, e.g., "CP" for pressure coefficient. 

Number of cases for which perturbation solutions are 
to be computed. 

Number of points (in addition to end points) to be 
held invariant in straining; note: 1 ^ NSELCT _< 6. 

Values of perturbation parameter in base, calibration, 
and perturbation solutions. 

Character string of length 80; identifies job and is 
printed as headline on first page of output. 
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XBASE (I ) , XCALB (I ) , XCHEK Cl ) . . . 


Arrays of surface coordinates in base, calibration, 
and comparison solutions. 

YBASE ( I ) , YCALB ( I ) , YCHEK ( I ) 

Arrays of dependent variables in base, calibration, 
and comparison solutions. 


A. 5 PREPARATION OF INPUT DATA 
A, 5.1 Description of Input 


Item 1 
Item 2 


Item 3 
Item 4 
Item 5 
Item 6 
Item 7 

Item 8 

Item 9 
Item 10 

Item 11 


One card, containing the parameters N, NCASE , LSPEC, 
LECHO, LUNIT, LCHEK, LPERT. 

One card, containing either 

(a) NSELCT, (LSELCT(I), I=1,NSELCT) 

(b) NSELCT, (LOCO (I), 1=1, NSELCT), 

(LOCI (I), 1=1, NSELCT) 

where (a) and (b) correspond to LSPEC = 0 and 
LSPEC = 1, respectively. 

One card, containing the character string TITLE. 

One card, containing the character string NAME. 

One card, containing the scaling parameters A and B. 
One card, containing MO (real) and QO. 

One set of K cards, where K = 1 + INT(N/8), containing 
data for x coordinate in base solution. 

One set of K cards, K as above, containing data for 
dependent variable in base solution. 

One card, containing Ml (real) and Ql . 

One set of K cards, K as above, containing data for 
X coordinate in calibration solution. 

One set of K cards, K as above, containing data for 
dependent variable in calibration solution. 
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Item 12 


One card, containing M2 (real) and Q2. 


Item 13 One set of K cards, K as above, containing data for 
X coordinate in comparison solution. This item is 
required only when LCHEK = 1, 

Item 14 One set of K cards, K as above, containing data for 

dependent variable in comparison solution. This item 
is required only when LCHEK = 1 . 

Note: Items 12-14 are required, in sequence, as many times as 

specified by NCASE. 
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A. 5.2 Format of Input Data 


w Item no. 1; 1 card 


Variable 

N 

NCASE 

LSPEC 

LECHO 

LUNIT 

LCHEK 

LPERT 


Card column 

5 

10 

15 

20 

25 

30 

35 


Format type 

I 

I 

I 

I 

I 

I 

I 


Item no. 2a 

(LSPEC = 0) ; 1 card 


/- LSELCT (NSELCT) 

Variable 

NSELCT 

LSELCT(l) 

LSELCT(2) 

— 





Card column 

5 

10 

15 

20 

25 

30 

35 


Format type 

I 

I 

I 

I 

I 





Item no. 2b (LSPEC =1): 1 card /- LOCO (NSELCT) LOCl(NSELCT) 


Variable 

NSELCT 

LOCO (1) 

— 

; — t 

LOCI (1) 

— 

— t 


Card column 

5 

10 

15 

20 

25 

30 

35 

1 

Format type 

I 

I 

I 

I 

I 

I 

I 

{ 


Item no. 3: 1 card 

Variable 
Card column 
Format type 

title 

10 

20 

30 

40 

50 

60 

70 

80 

A 


Item no. 4: 1 card 


Variable 
Card column 
Format type 


NAME 


Item no. 5 : 

1 card 


Variable 

A 

B 

7 

Card column 

10 

20 

( 

Format type 

1 

F 

1 































Item no. 6: 1 card 


Variable 
Card colvimn 
Format type 

MO 

QO 


1 

10 

20 


F 

F 



values per card 

Item no. 7 : 

K cards, K = 1 + INT(N/8), 8 

Variable 
Card column 
Format type 

XBASE (1) 

XBASE ( 2 ) 

XBASE ( 3 ) 

— 





1C 

20 

30 

40 

50 

60 

70 

80 

F 

F 

F 

F 

F 

F 

F 

F 

Item no, 8: 

K cards, K as above, 8 values per card 

Variable 
Card column 
Format type 

YBASE (1) 

YBASE (2) 

YBASE (3) 

— 





10 

20 

30 

40 

50 

60 

70 

80 

F 

F 

F 

F 

F 

F 

F 

F 

Item no. 9: 

1 card 

a lues per card 

Variable 
Card column 

Format type 
Item no. 10 

Ml 

. Q1 _ 


10 

20 

( 

F 

F 

J 

K cards, K as above, 8 v< 

Variable 
Card column 
Format type 

XCALB(i) 

XCALB(2) 

XCALB ( 3 ) 

— — 





10 

20 

30 

40 

50 

60 

70 

80 

F 

F 

F 

F 

F 

F 

F 

F 


Variable 
Card column 
Format type 


YCALB(l) 

YCALB ( 2 ) 

YCALB ( 3 ) 

; 





10 

20 

30 

40 

50 

60 

70 

80 

F 

F 

F 

F 

F 

F 

F 

F 

















































































U) 

cn 


Item no. 12 ; 1 card 


Variable 

M2 

Q2 







Card column 

10 

20 







Format type 

F 

F 

1 






Item no. 13 : 

K cards, 

K = 1 + INT (N/8) , 8 

values per 

card 




Variable 

XCHEK(l) 

XCHEK(2) 

XCHEK(3) 

— 





Card column 

10 

20 

30 

40 

50 

60 

70 

80 

Format type 

F 

F 

F 

F 

F 

F 

F 

P 

Item no. 14: 

K cards, 

K as above, 8 values per card 





Variable 

YCHEK(l) 

YCHEK{2) 

YCHEK(3) 

— 





Card column 

10 

20 

30 

40 

50 

60 

70 

80 

Format type 

F 

F 

F 

F 

F 

F 

F 

F 





































A. 6 DESCRIPTION OF OUTPUT 


The first output item consists of a banner page, and the card 
images of the input data, the latter only if LECHO = 1. 

The second item is a page headed by the job title, listing: 

1. the input parameters relevant to the actual calculation; 

2. the critical values of the dependent variable; 

3. the locations of the minimum, maximum, and critical 
points found by the program; 

4. the straining points selected; 

5. the invariant points. 

Results for unit straining of XBASE, and unit perturbation 
of the dependent variable are the third item output; this is done 
only if LUNIT = 1. 

The fourth item (repeated for each case computed) summarizes 
the results of the calculation. The Mach number, the value of the 
perturbation parameter, and the critical value of the dependent 
variable are printed first, followed by the locations of the 
minimum, maximum, and critical points in the perturbation solution 
and comparison solution (if any) . Then follows a table listing 
XBASE, YBASE, XCALB, YCALB, XPERT (the strained coordinate), and 
YPERT (the computed value of the dependent variable) . If LCHEK = 1, 
three additional columns list XCHEK, YCHEK, and YPERT (INT), the 
latter being interpolated values of YPERT (the computed solution) 
at the points given by XCHEK. This allows direct numerical com- 
parison of YPERT with YCHEK, since the values of XPERT and XCHEK 
do not coincide in general. 


A. 7 ERROR MESSAGES 

NUMBER OF CRITICAL POINTS IN 
BASE AND CALIBRATION SOLUTIONS 
ARE UNEQUAL - CALCULATION ENDED 

This message will be printed if critical points are specified 
in straining (LSPEC = 0) and the number of critical points in base 
and calibration solutions are unequal. The remedy is to avoid use 
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of critical points in straining, or to use base and calibration 
solutions having equal numbers of critical points. 


NUMBER OF CRITICAL POINTS 
SELECTED EXCEEDS NUMBER 
ACTUALLY LOCATED - CALCULATION 
ENDED 

This message will be printed if more critical points are 
specified in straining (LSPEC = 0) than the number located by the 
program. The remedy is to specify a number of points less than 
or equal to the actual number. 


ORDER OF SPECIFIED POINTS IN 
BASE AND CALIBRATION SOLUTIONS 
DOES NOT CORRESPOND 

This message will be printed if the fixed points specified 
(LSPEC = 0) occur in a different sequence in the base and calibra- 
tion solutions. The remedy is to use base and calibration solu- 
tions having the same qualitative features. 
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A. 8 SAMPLE CASE 


The sample case presented in this section provides results 
(6 perturbation calculations and comparisons with 'exact' 
nonlinear solutions) for a multiple-shock flow for which partial 
results were provided in figure 14 of the main text. The calcu- 
lation is for ahgle-of-attack perturbations of full potential 
flows past an isolated NACA 0012 airfoil at M^=0.80. The base 
and calibration angles-of-attack are aj 3 = 0.500° and 0.200°. 

Perturbation results are determined at a = 0.00°, 0.10°, 0.30°, 

0.40°, 0.60°, and 0.70° and are compared with previously-calcu- 
lated 'exact' nonlinear flows at those angles. 

The input data is tabulated in figure A.l, with item numbers 
corresponding to those indentified in Section A. 5.1 and A. 5. 2. 

The first card, item 0, indicates that there are 149 points 
(N=149) at which data will be input for the base, calibration, 
and comparison solutions; that there will be 5 cases (NCASE = 6) 
for which perturbation solutions are to be computed, that the 
invariant points will be located by the program (LSEPC=0), that 
the input card deck will not be printed (LECHO = 0), that the 
information regarding the unit perturbation will be printed 
(LUNIT = 1) , that there will be a comparison of the perturbation 
results with the exact solution (LCHEK=1), and that the plot 
output will denote an angle-of-attack perturbation CLPERT = 2 ) . 

The second card, item 2a, indicates that there will be three 
invariant points (NSELCT=3) in addition to the end points; and 
that those points will be (1) where the maximum occurs 
(LSELCT(l) = 2) i.e. the stagnation point, (2) the first critical 
point (LSELCT(2) =3) i.e. the 1st shock point found when moving 
forward on the bottom surface from the trailing edge, and (3) 
the second critical point (LSELCT(3) =4) i.e. the 2nd shock point 
The next card, item 3, contains the identifying title. On the 
next card, item 4, the 2 length character string indicates that 
the dependent variable for print output will be symbolized by a 
'CP' denoting pressure coefficient. Item 5 indicates that the 
coordinates of the data points to be read in will start at x=1.0 
on the upper surface (refer to descriptions in A.4). The next 
card, item 6, indicates that the base flow values of Mach number 
and perturbation parameter (angle-of-attack in this case) are 
MO = 0.80 and Q0=0.50, respectively. The following 19 cards, 

item 7, provide the 149 base flow values of the surface coordi- 
nates, while the next 19 cards, item 8, provide the 149 base flow 
values of the dependent variable (pressure coefficient) . Items 9 
10, and 11 indicate for the calibration flow the corresponding 
information given by the items 6,7, and 8 for the base flow. 

Items 12, 13, and 14, of which there are six sets corresponding 
to the 6 cases to be studied, provide analogous information as 
items 6,7, and 8, but now refer to the 'exact' nonlinear results. 
These, of course have been previously computed at the indicated 
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values of angle-of -attack (Q2) given in Item 12, and are 
included here for comparative purposed to enable assessment of 
the perturbation results. 

Figure A. 2 provides an abbreviated print output for the 
sample case, while figure A. 3 provides the plot output of the 
results for the six cases, and display the base (••••) ^ 

calibration ( ), perturbation (****), and 'exact' nonlinear 

( ) flow solutions. 
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ColTomn 

No. 


Item 

No. 


5 

H 

H 

to 

o 

r 1 

30 35 

t 1 

149j 

6 

0 0 

111 

ll 2 

3L 

2 

3 4 

r 


80 


SAMPLE CASE - ALPHA 

cpI 

'ERTURBATIpN 

FOR MULTI- 

SHOCK ISOLATED AIRFOIL FLOW 

. 






-1.0 1.0 
.6QQQ00 .500000 


.999495 .997758 

.9951 32 

.991597 

.987091 

.981656 

.975290 

.968014 

.959627 .950749 

.940799 

.930018 

.918425 

.906020 

.892863 

.878994 

.864393 .849140 

.633235 

.816758 

.799690 

.7621 1 0 

.764038 

.745536 

.726603 .707279 

.687646 

.667702 

.647510 

.6271 08 

.606518 

,685601 

.564955 .544062 

.523144 

.502239 

.481369 

.460595 

.439937 

.41941 7 

.399095 .378972 

.359110 

.339530 

.320253 

.301 322 

.282758 

.264582 

.246636 .229547 

.212731 

. 1 96413 

. 160616 

.165342 

.150655 

. 1 36558 

.123054 .110166 

.097919 

.066338 

.075427 

.065211 

.055678 

.046895 

.033773 .031422 

.024618 

.016960 

.013897 

.009586 

.006069 

.003379 

.001483 .000362 

. 000000 

.000362 

.001483 

.003379 

.006089 

.009586 

.015397 .018980 

.024818 

.031422 

.038773 

.046695 

.055678 

.065211 

.075427 .066338 

.097919 

,110166 

.123054 

. 1 36553 

.150655 

.165342 

.160616 .196413 

.212731 

.229547 

.246838 

.264582 

.282758 

.301 322 

.320253 .339530 

.359110 

.378972 

.399095 

.41 941 7 

.439937 

.460595 

.481369 .502239 

.523144 

.544062 

.564955 

.585301 

.606518 

.6271 08 

.647510 .667702 

.687646 

.707279 

.726603 

.745536 

.764038 

.7621 1 0 

.799690 .616758 

.833235 

.849140 

.664393 

.878994 

.692863 

.906020 

.918425 .930018 

.940799 

,950749 

.959827 

.968014 

.975290 

,981656 

.937091 .991597 

. 995132 

.997750 

.999495 




.494705 .494256 

.414440 

.401123 

.i70901 

: 33? 399 

.303686 

.274379 

.247635 .222349 

. 196566 

.174301 

.153539 

.133375 

. 110559 

.091755 

.072839 .053774 

.034598 

.016276 - 

.000869 

-.019323 

-.037707 

-.055942 

-.071935 -.090313 

-.109358 - 

.126665 - 

.146033 

-.164340 

-.183752 

-.202450 

-.221343 -.241733 

-.260654 - 

.278560 - 

.295525 

-.311698 

-.326422 

-. ;-j041 7 

-.403291 -.497623 

-.600835 - 

.668220 - 

. 700287 

-.71 1940 

-.712103 

-.707273 

-.699075 -.667197 

-.672616 - 

.656143 - 

,635506 

-.613957 

-.591617 

-.562885 

-.530617 -.496814 

-.460480 - 

.416394 - 

.372256 

-.319692 

-.268545 

-.205256 

-.131694 -.056739 

.037329 

.148457 

.285927 

.443489 

.629985 

.631 392 

1.019305 1.145394 

1.161699 1 

. 053470 

.849325 

.606187 

.371494 

.166265 

.005374 -.133558 

-.243041 - 

.335771 - 

.410394 

- .485396 

-.543403 

-.590331 

-.638336 -.679361 

-.716673 - 

.749187 - 

779261 

-.607499 

-.8341 31 

-.855738 

-.876312 -.896774 

-.914836 - 

.932147 - 

.947948 

-.961781 

-.974530 

-.986834 

-.997187 -1 .006749 

-1.016031 -1 

.022717 -1 

.029553 

-1 .034871 

-1 ,035277 

-1 .041615 

-1 .041697 -1 .039425 

-1.033151 -1 

.019966 - 

.993560 

-.947806 

-.857172 

-.470554 

-.047162 .060153 

.067065 

,059619 

.053938 

.049536 

.050777 

.055100 

.061931 .069486 

.079745 

.092114 

.105449 

. 1 1 9529 

,134216 

.153439 

.170576 .163732 

.210704 

.233000 

.256263 

.281298 

.309273 

,339810 

.374275 .403616 

.416177 

.495153 

.494705 




.eoboooj .2000001 







.996133 .994611 

.990123 

.984707 

.978365 

.971 116 

.962959 

.953914 

.944001 .933260 

.921710 

,909350 

.896242 

.662424 

.867877 

.852680 

.836834 .820417 

.80341 t 

.785896 

.767691 

.749457 

.730593 

.71 1341 

.691779 .671909 

.651791 

.631465 

610951 

.590310 

.569541 

.548725 

.527834 ,507056 

.486263 

.465566 

444934 

.424539 

.404292 

.384243 

.364455 ,344946 

.325740 

.306873 

.286381 

.270271 

.252590 

.235359 

.218602 .202540 

.186596 

.171371 

156730 

.142674 

.129207 

.116351 

.104130 .092568 

.081668 

071455 

061912 

.053105 

.044940 

,037523 

.030321 .024341 

.019557 

014959 

01 1060 

.007786 

.005125 

.002995 

.001399 .000360 

.000000 

00C360 

001399 

.002995 

.005125 

.007786 

.011060 .014959 

.019557 

024841 

030821 

.037323 

.044940 

.053105 

,061912 .071455 

.081668 

092568 

104130 

.1 16351 

.129207 

.142674 

.156730 .171371 

. 166596 

202340 

216602 

.235359 

.252590 

.270271 

.288381 .306678 

.325740 

344946 

364455 

.384243 

.404292 

.424539 

.444984 .465566 

.486263 

507056 

527884 

.548725 

.569541 

.590310 


Figure A.l- Card input for sample case 
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.610951 

.631465 

.651791 

.671909 

.691779 

.711341 

.730593 

.749457 

.767891 

.765896 

.803411 

.820417 

.636834 

.852660 

.867877 

.682424 

. 6962^2 

.909350 

.921710 

.933260 

.944001 

.953914 

.962959 

.971 116 

. 97 C 365 

.984707 

.990123 

.994611 

.998133 




.<+ 4916 *^ 

.448979 

.368724 

.356677 

.323122 

.288584 

.256968 

.226251 

.201234 

.174994 

.152770 

.131476 

.107828 

.088414 

.069020 

.049604 

.030206 

.011790 

-.005341 

-.023660 

-.041742 

-.059474 

-.074748 

-.092628 

-.109767 

-.125161 

-.141041 

-.156452 

-.170621 

-.181949 

-.190490 

-.196776 

- , 1 98872 

-.203294 

-.23571 9 

-.367492 

-.586069 

-.746810 

-.804612 

-.622557 

-.830010 

-.830430 

-.826029 

-.623224 

-.615167 

-.805731 

-.794610 

-.760953 

-.765371 

-.748500 

-.728200 

-.707375 

-.685933 

-.658961 

-.629807 

-.596386 

-.565587 

-.528007 

-.483403 

-.441165 

-.394459 

-.336198 

-.271604 

-.211304 

-.1 36609 

-.051358 

.052134 

.163559 

.296591 

.449001 

.621890 

.81 0551 

.994877 

1 . 1 29636 

1 .169160 

1 .096662 

.933160 

.727769 

.526208 

.346353 

.190742 

. 056809 

-.053975 

-.156994 

-.241464 

-.316023 

- . 57753 1 

-.443985 

-.500232 

-.542507 

-.563047 

-.626762 

-.662591 

-.693996 

-.723334 

-.7511 71 

-.777519 

-.793919 

-.819680 

-.639982 

-.857221 

-.673266 

-.886000 

-.901015 

-.912830 

-.923922 

-.932634 

-.940086 

-.946688 

-.949963 

-.952335 

-.951561 

-.946329 

-.937133 

-.917745 

-.885770 

-.760090 

-.450920 

-.135657 

-.055363 

-.048679 

-.054497 

-.059285 

-.05671 3 

-.056033 

-.048960 

-.039039 

-.030023 

-.017360 

-.003396 

.01 1570 

.025968 

.042113 

.059613 

.077446 

.095514 

.115310 

.136492 

. 156974 

. 1 78507 

.204143 

.230640 

.258908 

.289933 

.324329 

.357590 

.369390 

.449350 

.449164 




. 600000 J 

.OOOOGOI 







.999495 

. 997758 

.995132 

.991597 

.987091 

.981656 

.975290 

.968014 

.959827 

.950749 

.940799 

.930018 

.918425 

.906020 

.892663 

.878994 

.864593 

.849140 

.833235 

.816753 

.799690 

.7821 10 

.764038 

.745536 

.726603 

.707279 

.687646 

,667702 

.64751 0 

.627108 

.606518 

.585601 

.564955 

.544062 

.523144 

.502239 

,461369 

.460595 

.439937 

.419417 

.399095 

.378972 

.359110 

.339530 

.320253 

.301322 

.262756 

.264582 

.246838 

.229547 

.212731 

.196413 

.180616 

.165342 

.150655 

.136558 

.123054 

.110166 

.09791 9 

.086338 

.075427 

.06521 1 

.055678 

.046895 

.033773 

.031422 

.024818 

.013980 

.013897 

.009586 

.006089 

.003379 

.001433 

.000362 

.000000 

.000362 

.001483 

.003379 

.006089 

.009586 

.013897 

.018930 

.024818 

.031422 

.033773 

.046895 

.055678 

.065211 

.075427 

.086338 

.09791 9 

.1 10166 

. 1 23034 

. 1 36558 

.150655 

.165342 

.160616 

.196413 

.21 2731 

.229547 

.246658 

.264582 

.282758 

.301322 

.320253 

.339530 

.359110 

.378972 

.399095 

.41941 7 

.439937 

. 4-0595 

.481369 

.502239 

.523144 

.544062 

,564955 

.585801 

,606513 

.627108 

.647510 

.667702 

.687646 

.707279 

.726603 

.745536 

.764038 

.782110 

.799690 

.816758 

.833235 

.849140 

.664393 

.878994 

.892663 

.906020 

.918425 

.930013 

.940799 

-950749 

. 959327 

.966014 

.975290 

.981656 

.987091 

.991597 

.995132 

.997758 

,999495 




.436539 

.485539 

.406037 

. 392726 

.362350 

.326615 

.294714 

.265263 

.238623 

.213577 

.189264 

.164999 

. 1 4431 3 

.124278 

.101626 

.083098 

.064554 

.045974 

.027430 

.009934 

-.006143 

-.023236 

-.039876 

-.055862 

- . 068951 

-.084043 

-.097616 

-.108363 

-.118872 

-. 125296 

-.128739 

-.127030 

-.123011 

-.131038 

-.199763 

-.450296 

-.729836 

-.839648 

-.868610 

-.882942 

-.889736 

-.891484 

-.891 050 

-.836414 

-.880317 

-.872670 

-.862403 

-.851047 

-.838181 

-.822909 

-.605990 

-.787772 

-.766224 

-.744070 

-.721154 

-.692722 

-.662052 

-.629341 

-.593876 

- . 553475 

-.510373 

-.459592 

-.407829 

-.341769 

-.267177 

-.1 93761 

-.101301 

.003967 

.147792 

.309033 

.505174 

.724631 

.941542 

1 .107686 

1 .170402 

1 .107686 

.941542 

,724631 

.505174 

.309038 

.147792 

.008967 

-.101302 

-.193761 

-.267177 

-.341 769 

-.407829 

-.459592 

-.510373 

-.553475 

-.593876 

-.629341 

-.662053 

-.692722 

-.721154 

-.744070 

-.766224 

-.787772 

-.805991 

-.822989 

-.833182 

-.851047 
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.407381 

.427614 

.448046 

.468614 

.489297 

.510077 

.530690 

.551718 

.572520 

.593275 

.613903 

.634403 

.654716 

.674821 

.694678 

.714226 

.733466 

.75251 7 

.770739 

.788732 

.806256 

,823231 

.839636 

.855472 

.870659 

.8851 96 

.899005 

.912104 

.924455 

.935998 

.946732 

.956638 

.965677 

.973628 

.981072 

.937410 

.992622 

.997307 




.668728 

.668307 

.522300 

.441 943 

.369762 

,314595 

.273042 

.239321 

.209654 

.185547 

. 163233 

.139135 

.119414 

.099630 

.080298 

.060812 
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.042266 

.024934 

.006399 

-.011958 

-.099637 

-.1 18558 

-.136324 

-.155210 

-.252466 

-.273361 

-.295663 

-.316378 

-.437231 

-.471 741 

-.506264 

-.539430 

-.571404 

- .557257 

-.540964 

-.523283 

-.372779 

-.334533 

-.290600 

-.249155 

.06491 1 

.156050 

.250797 

.360544 

1.076692 

1 .162004 

1 .149582 

1 .027122 

! .049066 

-.077266 

-.187047 

-.260101 

-.619415 

-.666576 

-.705960 

-.744679 

-.833137 

-.906457 

-.925930 

-.94491 0 

- 1 .023092 

- 1 .035272 

- 1 .047226 

- 1 .057403 

- 1 .097319 

- 1.1 02275 

-1 . 1 06152 

- 1 . 1 11634 

- 1.11 7793 

- .112354 

- 1.098409 

- 1 .070169 

.154654 

.163569 

.160062 

.157352 

.160676 

.191367 

.207425 

.222816 

.382321 

.450996 

.528020 

.669150 


-.030082 

-.045955 

-.064578 

-.082762 

-.175483 

-.192160 

-.212669 

-.232627 

-.338889 

-.361301 

-.383949 

-.409232 

-.562679 

-.577161 

-.581608 

-.579150 

-.498949 

-.471497 

-.441655 

-.409089 

-.196666 

-.137268 

-.061064 

-.012045 

.482506 

.616924 

.770906 

.932586 

.829499 

.598496 

.360563 

.198798 

-.373953 

-.451 564 

-.515760 

-.564158 

-.777275 

-.808043 

-.634895 

-.859770 

-.963693 

-.979917 

-.995537 

- 1 .01 0093 

- 1.066934 

- 1 .076335 

- 1 .063463 

- 1 .090957 

- 1 . 114988 

-1 . 1 17839 

-1 . 119957 

- 1 .119530 

- 1 .009648 

-.890889 

-.417645 

.062612 

.156120 

.158521 

.164013 

.171487 

.240828 

.668728 

.265001 

.293844 

.331050 
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« « 

* PROGRAM PERTURB * 

w w 

* CALCULATES NONLINEAR SINGLE-PARAMETER * 

* » 

« CONTINUOUS OR DISCONTINUOUS * 

* PERTURBATION SOLUTIONS * 

* * 

* WHICH REPRESENT A CHANGE IN EITHER ft 

* w 

* GEOMETRY OR FLOW CONDITIONS » 

* * 

* BY EMPLOYING A STRAINED-COORDINATE PROCEDURE * 

* » 

« UTILIZING A UNIT PERTURBATION DETERMINED FROM » 

* N 

* TWO PREVIOUSLY CALCULATED * 

* 

* ’BASE' AND 'CALIBRATION' SOLUTIONS * 

* n 
DISPLACED FROM ONE ANOTHER BY SOME REASONABLE * 

* 

* CHANGE IN GEOMETRY OR FLOW CONDITION * 

» # 

* WRITTEN BY * 

* * 

^ JAMES P. ELLIOTT AND STEPHEN S, STAHARA * 

* * 

* # 

» NIELSEN ENGINEERING AND RESEARCH, INC. * 

» * 

« MOUNTAIN VIEW, CALIFORNIA « 




* SAMPLE CASE - ALPHA PERTURBATION FOR MULTI-SHOCK ISOLATED AIRFOIL FLOW # 


««<««< LIST OF INPUT PARAMETERS >>>>>>>>>> 


N = 



A = -1 .0 B = 1 

I .0 


BASE SOLUTION; 

MQ = 0.8000 

QO = 0.5000 

CALIBRATION SOLN: 

Ml = 0.8000 

Q1 = 0.2000 


Figure A. 2- Abbreviated print output for sanple case 
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««««« CRITICAL VALUES OF CP »»»»» 


BASE SOLUTION: CPCRIT = -0.^3^6 

CALIBRATION SOLN: CPCRIT = -0.43^f6 


««< LOCATIONS OF MIN., MAX., AND CRITICAL PTS. >»» 
DENOTES POINT ON LOWER SURFACE) 


BASE SOLUTION: 


MINIMUM AT X = 0.4606 
MAXIMUM AT X = 0.0000 
2 CRITICAL POINT(S): 
1ST AT X = 0.3924W 
2ND AT X = 0.6288 


(POINT #112) 

(POINT # 75) 

(AFTER POINT # 41 ) 
(AFTER POINT #120) 


CALIBRATION SOLN: 


MINIMUM AT X = 0.4043 
MA>IMUM AT X = 0.0000 
2 CRITICAL POINT(S) : 
1ST AT X = 0.4592^f 
2ND AT X = 0.5498 


( POINT #111) 

(POINT # 75) 

(AFTER POINT # 36) 
(AFTER POINT #118) 


««««« STRAINING POINTS SELECTED »»»»» 

NUMBER OF FIXED POINTS : 5 

FIXED POINTS SELECTED (IN ADDITION TO END POINTS) : 

POINT OF MAXIMUM CP 
CPCRIT (1ST OF 2) 

CPCRIT (2ND OF 2) 


««««« LOCATION OF FIXED POINTS »»»»» 


{* DENOTES POINT ON LOWER SURFACE) 

• BASE SOLUTION: 

XFIX( 1 ) = 1 .0000* 

XFIX(2) = 0.3924* 

XFIX(3) = 0.0000 
XFIX(4) = 0.6288 
XFIX(5) = 1 .0000 

CALIBRATION SOLN: 


XFIX( 1 ) 
XFIX(2) 
XFIX(3) 
XFIX(4) 
XFIX(5) 


1 . 0000 * 
0.4592* 
0.0000 
0.5498 
1.0000 
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» UNIT PERTURBATION OF CP » 
» C * 
« UNIT STRAINING OF XBASE * 



Figure A. 2- Continued 
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52 

0.1964 

53 

0.1806 

54 

0.1653 

55 

0.1507 

56 

0.1366 

57 

0.1231 

58 

0.1102 

59 

0.0979 

60 

0.0863 

61 

0.0754 

62 

0.0652 

63 

0.0557 

64 

0.0469 

65 

0.0388 

66 

0.0314 

67 

0.0248 

68 

0.0190 

69 

0.0139 

70 

0.0096 

71 

0.0061 

72 

0.0034 

73 

0.0015 

74 

0 . 00 Q 4 

75 

0.0000 

76 

0.0004 

77 

0.0015 

78 

0.0034 

79 

0.0061 

80 

0.0096 

81 

0.0139 

82 

0.0190 

83 

0.0248 

84 

0.0314 

65 

0.0388 

86 

0,0469 

87 

0.0557 

88 

0.0652 

89 

0 . 0/54 

90 

0.0863 

91 

0.0979 

92 

0.1102 

95 

0.1231 

94 

0.1366 

95 

0.1507 

96 

0.1653 

97 

0.1806 

98 

0.1964 

99 

0.21 27 

lOO 

0.2295 

101 

0.2468 

102 

0 . 2-46 

103 

0.2828 

104 

0.3013 

105 

0.3203 

106 

0.3395 

107 

0.5591 

108 

0.3790 

109 

0.3991 

110 

0.4194 

1 1 1 

0.4399 


0.2299 

0.3990 

0.2114 

0.4079 

0.1935 

0.4105 

0.1763 

0.4064 

0.1598 

0.4252 

0.1440 

0.4357 

0.1289 

0.4411 

0.1146 

0 . 4454 

0.1010 

0.4572 

0.0883 

0.4672 

0.0763 

0.4799 

0.0652 

0.4718 

0.0549 

0.4757 

0.0454 

0.4779 

0.0368 

0.4874 

0.0290 

0.4954 

0.0222 

0.4945 

0.0163 

0.51 34 

0.0112 

0.5078 

0.0071 

0.4605 

0.0040 

0.3529 

0.0017 

0.21 1 2 

0.0004 

0.0803 

0.0000 

- 0.0247 

0.0003 

- 0.1734 

0.0013 

- 0.3329 

0.0030 

- 0.4227 

0,0053 

- 0.4710 

0.0084 

- 0.4994 

0.0121 

- 0.4931 

0.0166 

- 0,5033 

0.0217 

- 0 . 491 1 

0.0275 

- 0.4720 

0.0339 

- 0.4489 

0.0410 

- 0.4620 

0.0487 

- 0.4517 

0.0570 

- 0.4046 

0.0659 

- 0.4024 

0.0755 

- 0.3963 

0.0856 

- 0.3821 

0.0963 

- 0.3693 

0.1076 

- 0.3594 

0-1194 

- 0.3552 

0.1317 

- 0,3520 

0.1446 

- 0.3368 

0. 1579 

- 0.3235 

0.1717 

- 0.3246 

0.1860 

- 0.3201 

0.2007 

- 0.3143 

0.2158 

- 0.3123 

0.2313 

- 0.3079 

0.2472 

- 0.3037 

0.2634 

- 0.3028 

0.2800 

- 0.2994 

0.2969 

- 0.2961 

0.3140 

- 0.2961 

0.3313 

- 0.2930 

0.3489 

- 0.2937 

0.3667 

- 0.2927 

0 . 3846 

- 0.2942 
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112 

0.4606 

0.4027 

- 0.2989 

113 

0.4314 

0.4209 

- 0.3000 

114 

0.5022 

0.4391 

- 0.3053 

115 

0.5231 

0.4574 

- 0.3079 

116 

0.5441 

0.4757 

- 0.3077 

117 

0.5650 

0.4939 

- 0.2921 

118 

0.5858 

0.5122 

- 0.2933 

119 

0.6065 

0.5303 

- 0.3833 

120 

0.6271 

0.5483 

- 0.0654 

121 

0.6475 

0.5724 

0.2949 

122 

0.6677 

0,5969 

0.3781 

123 

0,6876 

0.621 1 

0,3958 

124 

0.7073 

0.6449 

0,3910 

125 

0.7266 

0.6684 

0.3758 

126 

0 . 7455 

0.6913 

0.3521 

127 

0.7640 

0.7138 

0.3283 

126 

0.7821 

0.7357 

0.3056 

129 

0.7997 

0.7570 

0.2892 

130 

0.8168 

0.7777 

0.2640 

131 

0.8332 

0.7977 

0.2435 

132 

0.8491 

0.8170 

0.2301 

133 

0 . 6644 

0 . B 355 

0.2154 

134 

0.8790 

0.8532 

0.1976 

1 35 

0.8929 

0.8700 

0.1802 

1 36 

0.9060 

0.8660 

0.1773 

1 37 

0.9184 

0.9011 

0.1615 

1 38 

0.9300 

0.9151 

0.1423 

139 

0.9408 

0.9282 

0.1338 

140 

0.9507 

0.9403 

0.1259 

141 

0.9598 

0.9513 

0.1090 

142 

0.9680 

0.9612 

0.0929 

143 

0.9753 

0.9700 

0.0783 

144 

0.9817 

0.9777 

0.0614 

145 

0.9871 

0.9843 

0.0620 

146 

0.9916 

0.9898 

0.1164 

147 

0.9951 

0.9941 

- 0.0799 

148 

0 . 9978 

0.9973 

0.1531 

149 

0.9995 

0.9994 

0.1520 
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I 


« OUTPUT FOR CASE SI OF 6 * 

M2 = O.fiOOO 
Q2 = 0.0000 
CPCRIT = - 0 .^ 3^6 


««< LOCATIONS OF MIN., MAX., AND CRITICAL PTS. »»> 
(* DENOTES ^OINT ON LOWER SURFACE) 

PERTURB *T ION SOLN: 

MltJIMUM AT X = 0.^112* (POINT S ^^5 ) 

MAXIMUM AT X = 0.0000* (POINT # 75) 

2 CRITICAL POINT(S)! 

1ST AT X = 0.502^* (AFTER POINT # ^1) 

2ND AT X = 0.<i961 (AFTER POINT SI 20) 

COMPARISON SOLN! 

MINIMUM AT X = 0.3790* (POINT # ^2) 

MAXIMUM AT X = 0.0000 (POINT # 75) 

2 CRITICAL POINT(S): 


1ST 

AT 

X = 0.5035* 

(AFTER 

POINT 

# 35) 

2ND 

AT 

X = 0.5035 

(AFTER 

POINT 

#1 1A ) 


POINT 

XBASE 

CPBASE 

XCALB 

CPCALB 

XPERT 

CPPERT 

XCHEK 

CPCHEK 

CPPERK INT) 

1 

0.9995 

0.4947 

0.9981 

0.4492 

0.9996 

0.4189 

0.9995 

0.4865 

0.4189 

2 

0.9978 

0.4943 

0.9946 

0.4490 

0.9982 

0.4191 

0.9978 

0.4665 

0.4293 

3 

0.9951 

0.4144 

0.9901 

0.3687 

0.9960 

0.4721 

0.9951 

0.4060 

0.4556 


0.9916 

0.4011 

0.9847 

0.3567 

0.9931 

0.4186 

0.9916 

0.3927 

0.3946 

5 

0.9871 

0.3709 

0.9784 

0 . 3231 

0.9895 

0.3613 

0.9871 

0.3623 

0.3615 

6 

0.9817 

0.3354 

0.9711 

0.2684 

0.9650 

0.3618 

0.9817 

0.3266 

0.3433 

7 

0.9753 

0.3037 

0.9630 

0.2570 

0.9798 

0.3332 

0.9753 

0.2947 

0.3087 

6 

0.9680 

0.2744 

0.9539 

0 . 2283 

0.9739 

0,3010 

0.9680 

0.2653 

0.2750 

9 

0.9598 

0.2478 

0.9440 

0.2012 

0.9672 

0.2713 

0.9598 

0.2386 

0.2440 

10 

0.9507 

0.2228 

0.9333 

0. 1 750 

0.9598 

0.2438 

0.9507 

0.2136 

0.2155 

1 1 

0.9408 

0.1986 

0.9217 

0.1528 

0.9517 

0.2180 

0.9408 

0.1893 

0.1886 

12 

0.9300 

0.1743 

0.9093 

0.1315 

0.9428 

0.1936 

0.9300 

0.1650 

0.1638 

13 

0.9164 

0.1535 

0.8962 

0.1078 

0 . 9334 

0.1 705 

0.9184 

0.1443 

0.1416 

14 

0.9060 

0.1334 

0.8824 

0.0884 

0.9232 

0.1503 

0.9060 

0.1243 

0.1187 

15 

0.8929 

0.1106 

0.6679 

0.0690 

0.9125 

0.1313 

0.8929 

0.1016 

0.0972 

16 

0.8790 

0.0918 

0.8527 

0.0496 

0.9012 

0.1093 

0.8790 

0.0831 

0.0777 

17 

0.8644 

0.0728 

0.8368 

0.0302 

0.8893 

0 .091 9 

0.8644 

0.0646 

0.0581 

18 

0.8491 

0.0538 

0.8204 

0.0118 

0.8768 

0.0746 

0.8491 

0.0460 

0.0382 

19 

0.8332 

0.0346 

0.8034 

-0.0053 

0.8633 

0 . 0574 

0.8332 

0.0274 

0.0195 

20 

0.8168 

0.0163 

0 . 7859 

-0.0237 

0.8504 

0.0397 

0.8168 

0.0099 

0.0029 

21 

0.7997 

-0.0009 

0.7679 

-0.0417 

0.8564 

0.0227 

0.7997 

-0.0061 

-0.0151 

22 

0.7821 

-0.0193 

0.7495 

-0.0595 

0.8221 

0.0085 

0. 7821 

-0.0232 

-0.0338 

23 

0.7640 

-0.0377 

0.7306 

-0.0747 

0.8073 

-0 . 0072 

0.7640 

-0.0399 

-0.0496 

24 

0.7455 

-0.0559 

0.7113 

-0.0926 

0.7922 

-0.0228 

0.7455 

-0.0559 

-0.0645 

25 

0.7266 

-0.0719 

0.6918 

-0.1098 

0.7767 

-0.0396 

0.7266 

-0.0690 

-0.0818 
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26 

0.7073 

- 0.0908 

0.6719 

- 0.1252 

0.7609 

- 0.0520 

0.7073 

- 0.0840 

- 0.0961 

27 

0.6876 

- 0.1094 

0.6518 

- 0.1418 

0 . 7449 

- 0.0650 

0.6876 

- 0.0976 

- 0.1099 

28 

0.6677 

- 0.1267 

0.6315 

- 0.1564 

0.7286 

- 0.0803 

0.6677 

- 0.1084 

- 0.1225 

29 

0.6475 

- 0.1460 

0.6110 

- 0.1706 

0.7121 

- 0.0927 

0.6475 

- 0.1189 

- 0.1332 

30 

0.6271 

- 0.1643 

0.5903 

- 0.1819 

0.6955 

- 0.1043 

0.6271 

- 0.1253 

- 0.1385 

31 

0.6065 

- 0,1838 

0.5695 

- 0.1905 

0.6787 

- 0.1163 

0.6065 

- 0.1288 

- 0.1391 

32 

0.5858 

- 0.2024 

0.5487 

- 0.1968 

0.6617 

- 0.1259 

0.5858 

- 0.1270 

- 0.1356 

33 

0.5650 

- 0.2213 

0.5279 

- 0.1989 

0.6447 

- 0.1347 

0.5650 

- 0.1230 

- 0.1284 

34 

0.5441 

- 0.2417 

0.5071 

- 0.2033 

0.6277 

- 0.1385 

0.5441 

- 0.1311 

- 0.1 347 

35 

0.5231 

- 0.2608 

0.4863 

- 0.2357 

0.6106 

- 0.1394 

0.5231 

- 0.1998 

- 0.2159 

36 

0.5022 

- 0.2786 

0.4656 

- 0.3675 

0.5935 

- 0.1381 

0.5022 

- 0.4503 

- 0.4368 

37 

0.4814 

- 0.2955 

0.4450 

- 0.5861 

0.5765 

- 0.1327 

0.4814 

- 0.7298 

- 0.6896 

38 

0.4606 

- 0.3119 

0.4245 

- 0.7468 

0.5595 

- 0.1 263 

0.4606 

- 0.8396 

- 0.8651 

39 

0.4399 

- 0.3264 

0.4043 

- 0.8046 

0.5426 

- 0.1 355 

0,4399 

- 0.8686 

- 0.8992 

40 

0.4194 

- 0.3504 

0.3842 

- 0.8226 

0.5259 

- 0.1909 

0.4194 

- 0.8829 

- 0.9086 

41 

0.3991 

- 0.4033 

0.3645 

- 0.8300 

0.5093 

- 0.3436 

0.3991 

- 0.8897 

- 0.9096 

42 

0.3790 

- 0.4976 

0.3449 

- 0.8304 

0.4666 

- 0.6450 

0.3790 

- 0.8915 

- 0.9082 

43 

0.3591 

- 0.6009 

0.3257 

- 0.8280 

0.4611 

- 0-8644 

0.3591 

- 0.8910 

- 0 . 9056 

44 

0.3395 

- 0.6682 

0.3069 

- 0.8232 

0.4359 

- 0.9059 

0.3395 

- 0.8664 

- 0.9016 

45 

0.3203 

- 0.7003 

0.2884 

- 0.8152 

0.4112 

- 0.9100 

0.3203 

- 0.8803 

- 0.8942 

46 

0.3013 

- 0.71 19 

0.2703 

- 0.8057 

0.3869 

- 0.9091 

0.3013 

- 0.8727 

- 0.8360 

47 

0.2828 

- 0 . 71 21 

0.2526 

- 0.7946 

0.3630 

- 0.9064 

0.2828 

- 0.8624 

- 0.8765 

48 

0.2646 

- 0.7073 

0.2354 

- 0.7810 

0.3397 

- 0.9017 

0.2646 

- 0.8510 

- 0.8648 

49 

0.2468 

- 0.6991 

0.2186 

- 0.7654 

0.3169 

- 0.8929 

0.2468 

- 0.8382 

- 0.8514 

50 

0.2295 

- 0.6872 

0.2023 

- 0.7485 

0.2947 

- 0.8830 

0.2295 

- 0.8230 

- 0.8370 

51 

0.2127 

- 0.6726 

0.1866 

- 0.7282 

0.2731 

- 0.871 1 

0.2127 

- 0.8060 

- 0.8197 

52 

0.1964 

- 0.6561 

0. 1 714 

- 0.7074 

0.2522 

- 0.8556 

0.1964 

- 0.7878 

- 0.7995 

53 

0.1806 

- 0.6355 

0.1567 

- 0.6859 

0.2319 

- 0.8395 

0.1806 

- 0.7662 

- 0.7814 

54 

0.1653 

- 0.6140 

0.1427 

- 0.6590 

0.2123 

- 0.8192 

0,1653 

- 0.7441 

- 0.7601 

55 

0.1507 

- 0.5916 

0.1292 

- 0.6298 

0.1934 

- 0.7958 

0.1507 
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1^6 0.9<J16 

1^7 0.9951 

148 0.9978 

149 0.9995 


0.4036 0.9847 
0.4162 0.9901 
0.4952 0.9946 
0.4947 0.9981 
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0.3694 0.9934 
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0.4492 0.9993 
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0.3927 0.4144 
0.4060 0.4379 
0.4865 0.4186 
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* OUTPUT FOR CASE «2 OF 6 » 

M2 = 0.8000 
Q2 = 0.1000 
CPCRIT = -0.4346 


««< LOCATIONS OF MIN., MAX., AND CRITICAL PTS. >»» 
DENOTES POINT ON LOWER SURFACE) 

PERTURBATION SOLN: 

MINIMUM AT X = 0.3834 (POINT #112) 

MAXIMUM AT X = 0.0000* (POINT # 75) 

2 CRITICAL POINT(S): 

1ST AT X = 0.4805* (AFTER POINT t 41) 

2ND AT X = 0.5229 (AFTER POINT #120) 

COMPARISON SOLN! 

MINIMUM AT X = 0.4015 (POINT #110) 

MAXIMUM AT X = 0.0000 (POINT # 75) 

2 CRITICAL POINT! S)! 


1ST 

AT 

X = 

0.4816* 

(AFTER 

POINT 

» 36) 

2ND 

AT 

X = 

0.5258 

(AFTER 

POINT 

#116) 


POINT 

XBASE 

CPBASE 

XCALB 

CPCALB 

XPERT 

CPPERT 

XCHEK 

CPCHEK 

CPPERT! INT) 
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0.9995 

0.4947 
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0.4341 

0.9989 
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0.9765 
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26 

0.7073 

- 0.0908 

0.6719 

- 0.1252 

0.7502 

- 0.0598 

0.6895 

- 0.1068 

27 

0.6876 

- 0.1099 

0.6518 

“ 0.1918 

0.7335 

- 0.0738 

0.6696 

- 0.1209 

28 

0.6677 

- 0.1267 

0.6315 

“ 0.1569 

0.7169 

- 0.0896 

0.6999 

- 0.1396 

29 

0.6975 

- 0.1960 

0.6110 

“ 0.1706 

0.6992 

- 0.1039 

0.6291 

“ 0.1959 

30 

0.6271 

- 0.1693 

0.5903 

- 0.1819 

0.6818 

“ 0.1163 

0.6085 

- 0.1556 

31 

0.6065 

- 0.1838 

0.5695 

- 0.1905 

0.6692 

“ 0.1298 
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- 0.1610 

32 

0.5858 

“ 0.2029 

0.5967 
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0.6966 

- 0.1912 

0.5670 
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0.5279 
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0.6288 
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0.5962 
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C<<<<<<<<<<<<<<<<<<<<<<<< main program perturb >>»>»>>»>>»>»»»>>>| 

c I 

C CALCULATES CONTINUOUS OR DISCONTINUOUS NONLINEAR PERTURBATION I 

C SOLUTIONS WHICH REPRESENT A SINGLE-PARAMETER CHANGE IN EITHER I 

C GEOMETRY OR FLOW CONDITIONS BY EMPLOYING A STRAINED-COORDINATE I 

C PROCEDURE. THE METHOD UTILIZES A UNIT PERTURBATION, DETERMINEO I 

C FROM TWO PREVIOUSLY CALCULATED SOLUTIONS C’BASE' AND 'CALIBRATION' I 

C SOLUTIONS) OBTAINED FROM AN 'EXPENSIVE' COMPUTATIONAL PROCEDURE I 

C AND DISPLACED FROM ONE ANOTHER BY SOME REASONABLE CHANGE IN I 

C GEOMETRY OR FLOW VARIABLE, TO PREDICT NEW NONLINEAR SOLUTIONS OVERl 

C A RANGE OF PARAMETER VARIATION. I 

C I 

C THIS CURRENT VERSION OF THE PROCEDURE IS CONFIGURED TO PREDICT AND I 

C PLOT PRESSURE COEFFICIENTS ON A BUOE OR AIRFOIL SURFACE, AND CAN I 

C ACCOUNT FOR THE MOTION OFi I 

C (1) ONE OR MORE CRITICAL POINTS (SHOCK POINTS), I 

C (2) A STAGNATION POINT, I 

C (3) A MAXIMUM-SUCTION-PRESSURE POINT, I 

C OR SIMULTANEOUSLY FOR ANY COMBINATION OF THESE. I 

C i 

C THE PROGRAM IS ALSO CONFIGURED TO COMPARE THE PERTURBATION- I 

C PREDICTED SOLUTIONS WITH THE CORRESPONDING 'EXACT' SOLUTIONS 1 

r. OBTAINED BY EMPLOYING THE SAME 'EXPENSIVE' COMPUTATIONAL I 

C PROCEDURE USED TO DETERMINE THE BASE AND CALIBRATION SOLUTIONS. I 

C SEE THE SUBROUTINE INPUT FOR DETAILS. I 

C I 

C N = NO. OF POINTS IN SURFACE PRESSURE DISTRIBUTION - ASSUMED EQUALl 

C FOR BASE, CALIBRATION, AND PREDICTED DISTRIBUTIONS. I 

C NOTEi N <= 200. ( 


ONCOMING MACH NO. 
PARAMETER PERTURBED 


CALIBRATION 

Ml 

qi 


M2 = ONCOMING MACH NO. OF PREDICTED FLOW I 

Q2 = VALUE OF PERTURBED PARAMETER IN PREDICTED FLOW I 

COORDINATE STRAINING IS PIECEWISE LINEAR WITH END POINTS AND ONE I 
OB MORE USER-SELECTED INTERIOR POINTS HELD INVARIANT. I 

THE PROGRAM LOCATES MINIMUM, MAXIMUM, AND ALL CRITICAL POINTS 1 
(SHOCK POINTS) IN THE BASE AND CALIBRATION SOLUTIONS, AND STORES I 
THESE IN THE ARRAYS XLOCO AND XLOCt (IT IS ASSUMED THAT THE NUMBERl 
OF CRITICAL POINTS DOES NOT EXCEED FOUR) AS FOLLOWSi I 

BASE CALIBRATION I 


XLOCOd) = MINIMUM PT. 
XLOCO(2) = MAXIMUM PT. 
XLOCO(3) = CRITICAL PT. f1 


XLOCK 1 ) 
XLOCU2) 
XLOCI ( 3) 


MINIMUM PT. 
MAXIMUM PT. 
CRITICAL PT, •) 


XLOCO(6) = CRITICAL PT. 


CRITICAL PT. 


THE NUMBER OF POINTS SELECTED FROM THESE IS SPECIFIED BY N5ELCT. 
THE CORRESPONDING SUBSCRIPTS OF XLOCO AND XLOCI ARE SPECIFIED IN 
THE FIRST NSELCT ELEMENTS OF THE ARRAY LSECT, E.G. TO SELECT THE 
MAXIMUM POINT AND THE FIRST AND THIRD CRITICAL POINTS, ONE 
SPECIFIES: 


C LSELCT( 11=2 I 

C LSELCT(2) = 3 1 

C LSELCT( 31 = 5 I 

C I 

C PROVISION IS ALSO MADE TO ALLOW THE USER TO SPECIFY THE POINTS IN I 

C THE BASE AND CALIBRATION SOLUTIONS THAT HILL BE HELD INVARIANT, TQI 

C A MAXIMUM OF SIX. SEE THE SUBROUTINE INPUT FOR DETAILS. I 

C I 

DIMEMSTON XBASEI 200 ) ,XCALB( 200 I ,XCHEK( 200 ) , XPERT! 200 I ,XUN1T( 200 I I 
DIMENSION YBASE( 200 ) , YCALB( 200 ) , YPERT( 200 ) . YCHEK( 200 I , YINTP! 2001, I 
7. YPRTI(200),YUNIT( 200 I I 

DIMENSION XLOCO(6 ),XLOCI(6 ),XLOC2(6 I ,XLOC3( 6 ) ,XFIX0( 8 I ,XFIX1 ( 6 1 1 

DIMENSION LCRO(')),LCRI(‘>),LCR2(^),ICR3(^),LOCOI6),LOC1(6), I 

7 ISEG0(al,ISEqi(8l,LSELCT(6) I 

DIMENSION X0UT(8) I 

DIMENSION TITLE! 20) I 

IMTEGER»2 NAME I 

INTEGER HEAD0(5) /^HBASE.^H SOL,^iHUTIO,^HN! ,^H /, 1 

7 HEADl(B) /4HCALI ,AHBRAT,AHION .^HSOLN.^H: /, I 

7 HEAD2(5) /<iHPERT,^HURBA,'^HTION,^H SOL,<tHNi /, I 

7 HEAD3(5) /<,HC0MP.i(HARI5,AH0N S.^HOLN: ,:,H / I 

INTEGER ORDIAI /<iH1ST ,AH2ND ,4H3RD ,4H4TH / I 

LOGICAL*! FLAG! a ) I 

REAL M' ,H1 ,M2 I 

COMMON /COEFF/ C(7I,D(7) I 

COMMON /PARAM/ TITLE , LOCO , LOC I , LSE LCT ,M .NCASE , LSPEC > LECHO , LUHIT , 1 

7 LCHEK.LPERT, NSELCT, A, B, NAME I 

COMMON /PERT/ MO ,M1 , M2 , QO , QI , Q2 , YCRO , YCR 1 , YCR2 i 

COMMON /XY/ XBASE, XCALB, XPERT, XCH£K,YBASE,YCALB,YPERT,YCHEK I 

DATA LTEHM /O/, LCORR /O/ I 

C I 

C***********************************************************************! 
C I 

C USER-SUPPLIED STATEMENT FUtlCriOH YCRIT(Z) DETERMINES CRITICAL ) 

C VALUES OF FLOW VARIABLE YCRIT AS FUNCTION OF FLOW PARAMETER Z. I 

C IGRAO l+t OR -n IS THE USER-SUPPLIED ALGEBRAIC SIGN OF DYCflIT/DX I 

C USED IN LOCATING THE CRITICAL POINT. I 

C I 

C IN THE PRESENT VERSION OF THE CODE, YCRIT REPRESENTS THE FULL- I 

C POTENTIAL CRITICAL PRESSURE COEFFICIENT FOR AIR (GAMMA = 1.4), Z 1 

C IS THE FREE STREAM MACH NUMBER, AND IGRAO CORRESPONDS TO POSITIVE 1 

C PRESSURE GRADIENT (+1). I 

C I 

YCRIT(Z)=2.0»(((Z.O+0.4*Z»»2)/2.4I*«(1.4/0.4)-1.0)/().4»Z**2) I 

IGRAO =1 I 


C****************** **«»*»»****»***«***«* ************* ***!>*)(* ************ I 

c I 

C INPUT CONTROL, GEOMETRY, AND STRAINING PARAMETERS, AND DATA FOR I 

C BASE AND CALIBRATION SOLUTIONS. I 

C I 

CALL INPUT! 1 ) I 

C I 

C PRINT BANNER PAGE I 

C I 

CALL BANNER I 

C I 

C ECHO INPUT DECK IF LECHO .NE. 0, I 

C I 


IF (LECHO .NE. 01 CALL ECHINP 


c 


(MAIN12I 


XFIX1 1 1 + 1 )=XCALB( LOCI ( 1 ) ) 

IHAIN181 


DEL1=Qt-qO 

IMAIH122 


GO TO 10 

IMAINI 82 


tCRO = YCRIT(MO 1 

IMAINI23 


5 CONTINUE 

1 MAIM S3 


YCR1=YCRIT(m 1 

IMA1N124 


XFIX0(Y*1 ) = XLOCO(LSELCT(I>) 

IMAIN184 


WRITE 16,20001 TITLE 

1 MAIN) 25 


XFIXKI+I )=XL0C1(LSELCT(I)) 

IMAINI 85 


WRITE 16,20101 N,A.B,HEADO,nO,QO,HEADt ,M1 ,QI 

IMAINI26 


10 CONTINUE 

|HAlNt86 


WRITE <6,20201 NAME ,HEAOO. NAME, YCRO , HEAD) ,«AnE,TCRI 

1 MA1N127 


WRITE (6,2050) NFIX 

IMAIN187 

c 


IMAINI28 


IF (LSPEC .EQ. 0) GO TO 14 

IMAINI 86 

c.. 

...NORMALIZE X COORDINATES AND LOCATE MINIMUM, MAXIMUM, AND CRITICAL 

IMAIN129 


WRITE (6,2046) HEADO,HEADI 

IMAINI 89 

c 

POINTS FOR BASE AND CALIBRATION SOLUTIONS. 

1 MAIN! 30 


DO 13 1=1, NSELCT 

IMAIN190 

c 


IMAIN131 


WRITE (6.2047) LOCOd ) , LOCI ( I ) 

IMAINI 91 


CALL SCALE ( N , XBASE , ) , A ,B ) 

IMAINI32 


(3 CONTINUE 

IMAINI 92 


CALL LOCATE 1 H , XBASE , YBASE , YCRO . IGRAO , LMNO , LMXO ,NCR0 . LCRO ,XLOCO ) 

1 MAIN133 


GO TO IS 

1MAINI93 


CALL SCALE <N,XCALB,» ,A,B) 

1MAIH134 


Y4 CONTINUE 

1MAIN194 


CALL LOCATE ( N,XCALB,YCALB ,YCR1 , IGRAO , LMN1 , LMX1 ,NCR1 , LCRt >XL0C1) 

IMAIN135 


00 IS 1=1 .NSELCT 

IMAINI 95 


WRITE 16,20301 

IMAIN136 


IF (LSELCT(I) .EQ. 1) WRITE (6.2060) NAME 

IMAINI 96 


WRITE (6,2035) 

IMAIN137 


IF (LSElCT(I) .EQ. 2) WRITE (6,2070) NAME 

IMAINI 97 


CALL UPLOW ( A,B,XLOCO,6.NCRO+2,XOUT,FLAG) 

(MAIN 138 


IF (LSELCT1I) .LE. 2 ) GO TO 15 

IMAINI 98 


WRITE I6,20<i0) HEAOO,XOUT( 11 ,FLAGM ),LMNO,XOUT(2 1 ,FLAG( 2 >,LMXO 

IMAIN139 


LCORR=1 

IMAINI 99 


IF (NCRO .GT. Ot WRITE (6, 20^,5) NCRO, 

IMAIN140 


LPR=LSELCT( I )-2 

IHAIN200 


V. (0RD(I),X0UT(I«2),FLAG(I+2),LCR0(I»,I = 1 ,NCRO) 

1 MAIN! 41 


WRITE (6.2080) NAME >ORD( LPR ) .NCRO 

IMAIN201 


CALL UPLOW ( A,B,XLOCI ,6,NCRl*2,XOUT,FLAG) 



15 CONTINUE 

IMAIN202 


WRITE (6,20'<0l HEAOl ,XQUT( 1 l.FLAGC 11 ,LMNt ,XQUT12) >FLAG(21 .LMXt 

(MAIN143 


18 CONTINUE 

IMAIN2Q3 


IF (NCR1 .GT. 0) WRITE (6,2045) NCRI, 

IMAIN144 

C 


IMAIN204 


V. (0RDII),X0UT(It2),FLA6(It2),LCRI(II,I = 1,NCR1 ) 

IMAIN145 

C. 

ARRANGE SELECTED FIXED POINTS IN A MONOTONE SEQUENCE, 

IHAIN205 

c 


IMAIH146 

c 


IMAIN206 

c.. 

...CHECK FOR INVALID STRAINING SPECIFICATION IF LSPEC = 0. 

IMAIN147 


CALL SORT (NFIX.XFIXO.ISEOO) 

IMAIN207 

c 


inAIN148 


CALL SORT (NFIX.XFIXl ,ISEQ1) 

lMAIH20a 


IF 1 LSPEC .Eq. n GO TO 4 

IMAIH149 


WRITE (6.2090) 

IMAIN209 


ICOUNT=0 

IMAIN150 


WRITE (6,2035) 

IMAIN210 


DO 2 1=1 ,NSELCT 

IMAIN151 


CALL UPLOW 1 A,B,XFIXO,8,NFIX,XOUT,FLAG) 

IMAIN2II 


IF (LSELCT(I) .LE. 2) GO TO 2 

IMAINI52 


WRITE (6,2100) HEADO,(I,XOUT(I),FLAG(I),I=liNFIX) 

IMAIN212 


IC0UUT=IC0UNTt1 

1 MAIN! S3 


CALL UPLOW (A.B.XFIXl ,8,NFIX,X0UT,FLAG) 

IMAIN213 


IF (NCRO .NE. NCRI) LTERM=1 

1 MAIN 154 


WRITE (6,2100) HEAD! ,(I.X0UT(I).FLA6(I),I=1,NFIX) 

IMAIN214 


2 CONTINUE 

IMAIN155 

c 


IMAIN2I5 

c 


IMAIN156 

c. 

STOP EXECUTION IF ORDER OF OCCURRENCE OF CRITICAL POINTS IN BASE 

IMAIN2I6 

c. 

...STOP EXECUTION IF CRITICAL POINTS ARE TO BE USED IN STRAINING AND 

IMAIN157 

c 

AND CALIBRATION SOLUTIONS DOES HOT CORRESPOND. 

IMAIN217 

c 

NUMBER OF CRITICAL POINTS IN BASE AND CALIBRATION SOLUTIONS ARE 

IMAIN158 

c 


IHAIN2I8 

c 

UNEQUAL. 

IMAIN159 


IF (LSPEC .EQ. 1 ) GO TO 25 

IMAIN2I9 

c 


|MAim60 


DO 20 1=1 , NFIX 

IMAIH220 


IF ( LTERM .EO. 1 ) GO TO 900 

IMAINI61 


IF (ISEQOm .NE, ISEQKin GO TO 910 

IMAIN221 

c 


1MAIN162 


20 CONTINUE 

IMAIN222 

c. . 

...STOP EXECUTION IF NUMBER OF CRITICAL POINTS SELECTED EXCEEDS 

1 MAIM163 


25 CONTINUE 

IMAIH223 

c 

NUMBER ACTUALLY LOCATED. 

IMAIM164 

c 


IHAIN224 

c 


IMAINI65 

c. 

COMPUTE COEFFICIENTS IN UNIT STRAINING OF XBASE i 

IMAIN225 


IF IICOUHT .GT. NCRO) 60 TO 905 

IMA1N166 

c 


{MAIN226 

c 


IMAIN167 

c 

XSTR : C(I) ♦ DdHtXBASE, NSEG, 

IMAIN227 


4 CONTINUE 

IMAIM168 

c 


IMAIN228 

c 


IMAIN169 

c 

WHERE NSEG IS THE NUMBER OF LINEAR SEGMENTS. 

IMAIN229 

c. 

...LOAD SELECTED STRAINING POINTS INTO FIXED-POINT ARRAYS FOR BASE 

1 MAIM 70 

c 


IMAIN230 

c 

AND CALIBRATION SOLUTIONS. 

1MAIN171 


NSEG=NFIX-1 

IMA1N231 

c 


IMAINI 72 


DO 30 1=1 .NSEG 

IHAIH232 


NFIX^NSELCI+2 

IMAINI 73 


CMUM=XFIX1 d )«XFIX0d*1 )-XFIXld+l )»XFIX0d) 

IMAIN233 


XFIX0( ( ) = 0.0 

IMAINI 74 


DMUM = XFIXI( I+l )-XFIX1d) 

IMAIN234 


XFIXK I 1 = 0.0 

IMAINI 75 


DENOM = XFIX0dtl )-XFIX0d) 

IMAIN235 


XnX0IUFIX|:| .0 

IMAINI 76 


C( I ) = CNUM/DENOM 

IMAIN236 


>:fix((( fix)=i .0 

IMAIN177 


D( I l=DNUM/DENOM 

IMAIN237 


DO 10 : 1 .NSELCT 

ImAINI 78 


30 CONTINUE 

IMAIN238 


IF 1 LSt'EC .EQ. 0 ) GO TO 5 

iMAim 79 

c 


IMAIN239 


XF IXOl : ♦ 1 ) = XBASEl LOCO! I ) ) 

iMAIHieO 

c. 

DETERMINE UNIT STRAINING OF XBASE 

tMAlH240 


cn 

Ln 



cr» 

o\ 







IMAIN24I 

C 


IMAIN301 



CALL STRAIN ( N.NSEG.XFIXO .XBASE , 1 . 0 .XUNIT ) 



IMAIM242 

C. 

LOCATE MINIMUM, MAXIMUM, AND CRITICAL POINTS IN PERTURBATION- 

IMAIN302 






IMAIN243 

C 

SOLUTION. 

IMAIN303 

c. 


INTERPOLATE CALIBRATION SOLUTION TO BASE FLOW POINTS CORRESPONDING! t1AIN2^4 

c 


|MA1N304 

c 


TO UNIT STRAINING. 



IMAIN245 


CALL SCALE ( N.XPERT, 2 , A ,B ) 

IMAIN505 

r 





IMAIH246 


CALL SCALE ( N.XPERT, 1 , A ,B ) 

IMA1N306 



CALL INTERP (N.XCALB.YCALB.XUNIT.YINTP) 



IMAIN247 


CALL LOCATE ( N.XPERT, YPERT, YCR2 , IGRAD , LHN2 , LMX2 ,NCR2 , LCR2 ,XLOC2 ) 

IMAIN307 

r 





|MAIN24a 


WRITE (6,2130) ICASE .NCASE ,H2 ,02 .NAME . YCR2 

IMAIN308 

c. 


CORRECT VALUES ON EITHER SIDE OF CRITICAL POINTS, IF THESE ARE 

IMAIN249 


WRITE (6,2030) 

IMAIN309 

c 


USED IN STRAINING. 



IMAIN250 


WRITE (6.2035) 

IHAIN310 






IMAIN251 


CALL UPLOW ( A,e,XLOC2.6,NCR242,XOUT,FLAG) 

IMAIN3) 1 



IF (LCORR -EQ. 0) GO TO 36 



IMA1N252 


WRITE (6,2040 ) HEAD2 ,XOUT( 1 ) .FLAGM ) , LMN2 ,XOUT( 2 ) ,FLAG( 2 ) , LMX2 

IHAIN312 



DO 35 1=1 .NCRI 



IMAIN253 


IF (NCR2 .GT. 0) WRITE (6,2045) NCR2, 

IMAIN313 



YINTPI LCR01 1 1 ) = YCALB( LCRU I )) 



IMAIN254 


•/. (0RD(I),X0UT(I42),FLAG(I42),LCR2(I),I = 1 ,NCR2 ) 

IMAIN3I4 



YINTP1 LCROlIl + 1 )=YCALB(LCR1(I) + 1 1 



IMAIN255 


CALL SCALE (N.XBASE ,2 , A ,B ) 

IMAIN3I5 


35 

CONTINUE 



IMAIN256 

c 


IMAIN3I6 


16 

CONTINUE 



IMAIN257 

c. 

....IF LCHEK .NE. 0 READ IN DATA FOR COMPARISON SOLUTION AND LOCATE 

IMAIN317 

c 





iMAINESa 

c 

MINIMUM, MAXIMUM, AND CRITICAL POINTS. 

IMAIH316 

c. 


RESTORE PHYSICAL X IN CALIBRATION SOLUTION SINCE IT HILL 

NOT 

BE 

IMAIN259 

c 


IMAIN319 

c 


FURTHER USED. 



1 MAIN260 


IF (LCHEK .EQ. 0 ) GO TO 90 

IMAIN320 

c 





IHAIN261 


CALL IHPUT<3) 

IMAIN321 



CALL SCALE (N,XCALB,2,A,8) 



IMAIN262 


CALL SCALE (N.XCHEK.I ,A,B) 

IMAIN322 

r 





IMAIN263 


CALL LOCATE ( N .XCHEK , YCHEK , YCR3, IGRAD , LMN3, LMX3 ,NCR3 , LCR3.XLOC3 ) 

IMA1N323 

c. 


DETERHrNE THE UNIT PERTURBATION. 



IMAIN264 


CALL UPLOW (A,B.XLOC3,6,NCR3*2,XOUT,FLAG) 

IMAIN324 






IMA1N265 


WRITE (6,2040) HEAD3,XOUT( 1 ) ,FLAG( 1 ) , LHN3,XOUT( 2 ) ,FLAG( 2 ) , LMX3 

IMA1N325 



DO 40 1=1 .N 



IMAIN266 


IF (NCR3 .GT. 0) WRITE (6,2045) NCR3, 

IMAIN326 


<iO 

YUNIK I ) = ( YINTPC 1 1-YBASEI I ) 1/DEL1 



IMA1N267 


y. (0RD(I),X0UT(I42),FLAG(I42),LCR3(I).I = 1 ,NCR3) 

IMAIN327 

r 





|riAlN26S 


CALL INTERP ( N.XPERT, YPERT, XCHEK, YPRTI ) 

IMAIN328 

c. 


.PRINT UNIT PERTURBATION AND UNIT STRAINING IF LUNIT .NE. 

0. 


IHAIN269 


CALL SCALE ( N.XPERT ,2 , A.B ) 

IMAIN329 

f 





IMAIN270 


CALL SCALE (N.XCHEK ,2 , A ,6 ) 

IHAIH330 



IF (LUNIT .EQ. 0) GO TO 50 



IHAIN271 


WRITE (6,2135) NAME .NAME .NAME .NAME .NAME 

1MAIN331 



CALL SCALE ( N.XBASE ,2 ,A>B 1 



IMAIH272 


WRITE (6.2140) (I,XBASE(I),YBASE(I).XCALB(I),YCALB(I), 

IHAIN332 



CALL SCALE ( N.XUNIT ,2 ,A>B ) 



IMAIH273 


y. XPERT ( I ) , YPERT( I ) , XCHEK( I ) , YCH£K( I ) , YPRTK I ) , 

1MAIN333 



WRITE (6,2110) NANE.NANE 



IMAIN274 


y. 1 = 1 , N) 

IHAIN334 



WRITE (6,2120) ( I ,XBASE( I ( ,XUNIT( I ) ,YUHIT( I ) ,1 = 1 ,N ) 



IMAIN275 


GO TO 100 

IMAIN335 



CALL SCALE ( N,XBASE . 1 . A.B ) 



IMAIN276 


90 CONTINUE 

IMAIN336 


50 

CONTINUE 



IMAIN277 

c 


IMAIN337 

c 





|MAIH27a 


CALL SCALE ( N.XPERT , 2 , A ,B ) 

IMAIN338 

c. 


■CONSTRUCT PERTURBATION SOLUTIONS FOR TEST CASES (AND COMPARE 

WITH 

1MAIN279 


WRITE (6,2145) NAME .NAME .NAME 

IMAIN339 

c 


EXACT SOLUTION, IF AVAILABLE). 



IMAIN280 


WRITE (6,2150) II,XBASE(I),YBASEII),XCALB(I),YCALB(I), 

IMA1N340 






IHA1N2S1 


y. XPERT) I),YPERT(I),I = I ,N) 

IHAIN341 



DO 200 ICASE=1 ,NCASE 



|MAIN2e2 


100 CONTINUE 

IMAIN342 



CALL INPUT! 2) 



|MAIN263 

C 


IHAIN343 



DEL2=Q2-qO 



lMAIN2a4 

c 

IF LCHEK .NE. 0 PLOT PERTURBATION AND COMPARISON SOLUTIONS. 

IMAIN344 



OEL21=DEL2/DELl 



IMAIN2B5 

c 


IHA1H345 



YCR2=YCRIT(M2) 



iMAINEea 


IF (LCHEK .EQ. 11 CALL PLOT (N.LPERT) 

IMAIN346 



YCR3=YCR2 



IMAIM287 


CALL SCALE ( N .XBASE , 1 , A.B ) 

IHAIH347 

c 





IMAIN288 


200 CONTINUE 

|MAIN34a 

c. 


.DETERMINE STRAINED COORDINATE FOR GIVEN PERTURBATION. 



IMAIN289 


GO TO 999 

IMAIN349 

c 





IMAIN290 

c 


{MAIN 550 



CALL STRAIN ( N,NSEG,XFIXO , XBASE ,OEL2) .XPERT ) 



IMAIN29I 

c 

ABNORMAL TERMINATION OF COMPUTATION. 

IMAIN351 

c 





IMA1H292 

c 


IMA1N352 

c. 


.DETERMINE PERTURBATION SOLUTION. 



IMAIN293 


900 WRITE (6,9000) 

IMA1N353 

c 





IMAIN294 


GO TO 9 99 

IMAIN354 



DO 60 1=1 ,N 



{MAINE 95 


905 WRITE (6,9050) 

IMAIN355 


60 

YPERT1 1 ) = YBASE( I )+DEL2»YUNIT( I ) 



IMAIN296 


GO TO 99 

IMAIN356 

c 





IMAIN297 


910 WRITE 16,9100) 

II1AIN3S7 

c. 


.ADJUST VALUES NEAR CRITICAL POINT FOR MONOTONE BEHAVIOR. 



IHAIN298 

c 


|MAIN35a 

c 





IMAIN299 


999 WRITE (6,9500) 

IMAIN369 



IF (LCORR .EQ. 1) CALL MONO ( NCRO, LCRO .XPERT, YPERT ) 



IMAIN300 


STOP 

IMAIN360 



c 

C I/O FORMAT STATEMENTS FOLLOH. 

C 

tOOO FORMAT (tHI il3Z(1H«>/ 

•/. 1X>1H«>ZSX,20A4i25X,1H«/ 

•/. IX,132( 

2010 FORMAT ( IH « lOUHO , 1X.24HLIST OF INPUT PARAMETERS, IX, 1 0MH> )// 
y, 6X,3HN =,tX,I3// 

y. 6X,3HA =.1X,F<*.1,<»X,3HB =,1X,F4.t// 

y. 6X,5A4,1X,AHM0 =,1X,F6.A,AX,AHqO =,1X,F6.4//’ 

y. 6X,5A<i,1X,4HM1 =,1X,F6.4,4X,4HQ1 =,1X,F6.4///» 

2020 FORMAT MH , 1 0MH< i , IX, 18KCRIT1CAL VALUES OF, IX, A2 . IX, 1 0(1H> )// 
y. 2(6X,5A4,IX,A2,6HCRIT =,IX,F7.4//I/) 

2030 FORMAT tIH ,5UH<).IX.37HL0CATI0NS OF MIN., MAX., AND CRITICAL, 
y. IX,4HPTS.,IX,5IIH>n 
2035 FORMAT (IH ,2X,34H(» DENOTES POINT ON LONER SURFACE)) 

2040 FORMAT (/6X,5A4// 

y. IIX,I4HMINIMUM AT X =, 1X,F6 .4,AI ,3X,8HI POINT t,I3,IH)/ 

y. IIX,|4HHAXIMUM AT X =,IX,F6.4.AI ,3X,8HIP01NT t,I3,IH)) 

2045 FORMAT (IH ,tOX,II ,IX,I8HCR1TICAL POINTISH/ 

X ( 15X,A4,6HAT X =,IX,F6.4,AI ,3X, 

y. I4H( AFTER POINT 0,13, IH)) I 

2046 FORMAT (IH , tOX,2( 5A4 I/I 

2047 FORMAT (IH , I4X,2HX( ,13, IH ), I5X,2HX( ,13, IH 1 1 

2050 FORMAT I///IX, 101 IH< I , IX,25HSTRAIN1NG POINTS SELECTED, IX, 

X 10MH>),// 

X 6X,24KMUMBER OF FIXED POINTS i,1X,I1// 

X 6X,51HFIXED POINTS SELECTED (IH ADDITION TO END POINTS) i 

X /) 

2060 FORMAT (IH ,10X.16HPOINT OF MINIMUM, I X,A2) 

2070 FORMAT ( IH ,IOX,I6HPOINT OF MAXIMUM, IX, A2 I 
2080 FORMAT ( IH ,I0X,A2,6HCRIT ( ,A4,3H0F ,II,IHI) 

2090 FORMAT (///IX, t0( IH< ) , IX,24HLOCATION OF FIXED POINTS, IX, IO( IH> I ) 
2100 FORMAT l/6X,5A4// 

X MH ,I0X,5KXFIX(,II,3H) s,IX,F6.4,AI 1 1 

tlfO FORMAT MHI,27MH)))/ 

X 1X,IH«,IX,20HUNIT PERTURBATION OF , IX, A2, IX, IH*/ 

X tX,1H«,t2X,IHC,l2X,IH»/ 

X tX,1H»,IX,23HUNIT STRAINING OF XBASE , IX, IH»/ 

X IX,27MH4)/// 

X I X , 5HPOINT , 4X , 5HXB ASE , 4X , 8HXSTRUNIT , 3X , A2 , 4HUNIT/ ) 

2120 FORMAT (IH ,1X,I3,IX,3F10.4I 
2130 FORMAT ( 1HI ,27< 1H«)/ 

X IX,I9HM OUTPUT FOR CASE 0,11 >4H OF ,II,2H «/ 

X IX,27(IH»I// 

X 6X,4HM2 =,1X,F6.4// 

X 6X,4HQ2 =,IX,F6.4// 

X 6X,A2,6HCRIT =,IX,F7.4///I 

2135 FORMAT (///IX,5HPOINT,4X,5HXBASE,5X,A2,4HBASE, 

X 4X,5HXCALB,5X,A2,4HCALB, 

7. 4X,5HXPERT,5X,A2,4HPERT, 

X 4X,5HXCHEK,5X,A2,4HCHEK, 

X 2X,A2,9HPERT(INTI/I 

2140 FORMAT ( IH , IX,I3,IX,9F10.4I 
2145 FORMAT (///IX,5HPOINT,4X,5HXBASE,5X,A2,4HBASE, 

X 4X,SHXCALB,5X,A2,4HCALB, 

X 4X,5HXPERT,5X,A2,4HPERT/I 

2150 FORMAT I IH ,|X,I3,IX,6F10.4) 

9000 FORMAT (///IX,26HNUHBER OF CRITICAL POINTS IN/ 

X IX,30HBASE AND CALIBRATION SOLUTIONS/ 


•vl 


X IX,3IHARE UNEQUAL - CALCULATION ENDED) 

9050 FORMAT (///IX.25HNUMBER OF CRITICAL POINTS/ 

X 1X.23HSELECTED EXCEEDS lAJHBEfi/ 

X IX.30HACTUALLT LOCATED - CALCULATION/ 

X IX,SHENDED) 

9100 FORMAT (///1X,28HOHDER OF SPECIFIED POINTS IN/ 

X 1X,30HBASE AND CALIBRATION SOLUTIONS/ 

X IX,39HD0ES NOT CORRESPOND - CALCULATION ENDED) 

9500 FORMAT MHI > 

END 

SUBROUTINE INPUT (ICALL) 

C 

C ALL INPUT FOR PROGRAM PERTURB IS READ Bt THIS SUBROUTINE, AND IS 
C REQUIRED IN THE FOLLOUING ORDER (FOR DETAILS. REFER TO 
C ACCOMPANYING MANUAL). 


CARO tl (1615 1 
C 


NUMBER OF DATA POINTS IN BASE AND CALIBRATION I 

SOLUTIONS. I 

NUMBER OF CASES FOR HHICH PERTURBATION SOLUTIONS ARE I 
TO BE COMPUTED. I 

CONTROLS HOM INVARIANT POINTS IN STRAINING ARE 1 

SPECIFIED (SEE CARD *2). I 

LSPEC : 0 ... INVARIANT POINTS SELECTED FROM AMONG I 
THOSE LOCATED BY THE PROGRAM I 

LSPEC = 1 . . . INVARIANT POINTS PRESELECTED BY USER 1 

CONTROLS HHETHER OR NOT INPUT DECK IS PRINTED. I 

LECHO = 0 ... NO OUTPUT I 

LECHO = ( ... OUTPUT I 

CONTROLS HHETHER OR HOT UNIT COORDINATE STRAINING I 
AND UNIT PERTURBATION ARE PRINTED. I 

LUNIT = 0 ... NO OUTPUT I 

LUNIT = 1 ... OUTPUT I 

SPECIFIES HHETHER OR NOT PERTURBATION SOLUTION IS TOl 
BE CHECKED UITH AND PLOTTED AGAINST AH EXACT I 

COMPARISON SOLUTION. I 

LCHEK = 0 ... NO COMPARISON, NO PLOT I 

LCHEK = I ... COMPARISON, PLOT I 

SPECIFIES TYPE OF PERTURBATION. OPERATIONAL ONLY I 

WHEN LCHEK = 1 AND AFFECTS ONLY OUTPUT FROM PLOT I 

SUBROUTINE. I 

LPERT : I ... THICKNESS-RATIO PERTURBATION I 

LPERT = 2 ... ANGLE-OF-ATTACK PERTURBATION I 

LPERT : 3 ... MACH-NUM&ER PERTURBATION I 


C*»»)t CARD »2 (1615) «»in*)«)()n*»HHHH«)nnm)*)HHH*«)*)t)i) 
C 

C REQUIRED INPUT DEPENDS ON VALUE OF LSPEC > 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


xzmzzxzzxzxzKzzzzxzs: input fop lspec = o zxxzxxxxxzxzxxzzzxxxxxx 

•/. 

NUMBER OF POINTS (IN ADDITION TO END POINTS) TO BE 
HELD INVARIANT IN STRAINING, NOTEi 1 <= NSELCT <= 6 


NSELCT 


LSELCT(I) 


ARRAY OF LENGTH 6 OF WHICH NSELCT ELEMENTS ARE READ 
IN. SPECIFIES NATURE OF POINTS TO BE HELD INVARIANT 
ACCORDING TO THE CODE> 


MINIMUM PT. HELD INVARIANT 
MAXIMUM PT. HELD INVARIANT 
1ST CRITICAL PT. HELD INVARIANT 
2ND CRITICAL PT. HELD INVARIANT 
3RD CRITICAL PT. HELD INVARIANT 
4TH CRITICAL PT. HELD INVARIANT 


NOTE THAT THE CODE NUMBERS CAN BE ASSIGNED IN ANY 
ORDER. E.G. 


V. 

y. 

V. 

X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXXXXXXXXXXXX INPUT FOR LSPEC = 1 XXXXXXXXXXXXZXZXXXXXXXX 
X 

NUMBER OF POINTS (IN ADDITION TO END POINTS) TO BE 
HELD INVARIANT IN STRAINING. NOTE* 1 <= NSELCT <= 6. 


IS E(WIVALENT TO 


LSELCTd ) 
LSELCT(2) 
LSELCT(3) 


LSELCTt 1 ) 
LSELCT(2) 
LSELCT(3) 


BOTH CORRESPONDING TO NSELCT = 3 WITH THE MINIMUM, 
AND FIRST AND SECOND CRITICAL POINTS HELD INVARIANT. 


NSELCT 


ARRAY OF LENGTH 6 OF WHICH NSELCT ELEMENTS ARE READ 
IN. SPECIFIES SUBSRIPTS OF THOSE BASE FLOW POINTS 
WHICH ARE TO BE HELD INVARIANT. 

ARRAY OF LENGTH 6 OF WHICH NSELCT ELEMENTS ARE READ 
IN. SPECIFIES SUBSCRIPTS OF THOSE CALIBRATION FLOW 
POINTS WHICH ARE TO BE HELD INVARIANT. 


X 
X 
X 

X LOCO(I) 

X 
X 
X 

X LOCKI) 

X 
X 
X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxzzxxx 

CARD >3 (20A4) 

TI'^LE IDENTIFIES JOB - PRINTED AS HEADLINE ON FIRST PAGE 
OF OUTPUT. 

CARD (A2) 

NAME CHARACTER STRING OF LENGTH 2 WHICH SYMBOLIZES 

DEPENDENT VARIABLE, E.G. 'CP' FOR PRESSURE 


lINPUOSI 
I INPU0S2 
IINPU053 
I INPU054 
IINPU055 
IINPU056 
I INPU057 
I INPU058 
I INPU059 
I INPU060 
IINPU06I 
I INPU062 
I INPU063 
I INPU064 
I INPU065 
I INPU066 
I INPU067 
I 1NPU068 
I 1NPU069 
I INPU070 
IINPU071 
IINPU072 
IINPU073 
|INPU07<» 
IINPU075 
I INPU076 
I INPU077 
I INPU07B 
I INPU079 
I INPUOaO 
lINPUOSt 
I INPU002 
I INPU083 
INPU084 
INPUOaS 
iNPUoaa 
INPU0B7 
iNPUoas 
INPU0a9 
INPU090 
INPU091 
INPU092 
INPU093 
INPU099 
INPU095 
INPU096 
INPU097 
INPU098 
INPU099 
INPU100 
INPU101 
INPU102 
INPUI03 
1NPU104 
INPU105 
INPUI06 
INPU107 
INPUlOa 
INPU109 
INPU110 


C COEFFICIENT. llNPUlIl 

C IINPUI12 

CtiMMN CARD ts (aFIO.61 

c IINPU1I4 

C A SCALING PARAMETER (A = -X(1), WHERE X( II IS FIRST lINPUtlS 

C DATA POINT ON LOWER SURFACE ... SEE MANUAL). IINPUII 6 

C IINPUM7 

C B SCALING PARAMETER (D = X(N), WHERE X(NI IS LAST DATAIINPUII 8 

C POINT ON UPPER SURFACE ... SEE MANUAL). |INPUt19 

C IINPUI20 

C**4*iN CARO 86 (8F10.6) INPUt 21 

C IINPUI 22 

C MO ONCOMING MACH NUMBER IN BASE SOLUTION. IINPU123 

C IINPU124 

C QO VALUE OF PERTURBATION PARAMETER IN BASE SOLUTION. IINPU125 

C IINPU126 

ONE SET OF K CARDS (8F10.6), WHERE K = 1 ♦ INT(N/8) *****«**»H*tHHt*n*iii INPU1 27 
C |INPUt28 

C XBASE(I), 1=1, N ... IINPU129 

C X COORDINATE IN BASE SOLUTION. IINPU130 

C IINPU131 

ONE SET OF K CARDS (6F10.6), K AS ABOVE INPU1 32 

C IINPU133 

C YBASE(I), 1=1, N ... IINPU134 

C DEPENDENT VARIABLE IN BASE SOLUTION. IINPU135 

C IIHPU136 

c«««« NEXT CARD (8FI0.6) INPUl 37 

C IINPU138 

C Ml ONCOMING MACH NUMBER IN CALIBRATION SOLUTION. I1NPU139 

C IINPU140 

C (J1 VALUE OF PERTURBATION PARAMETER IN CALIBRATION IlNPUtAl 

C SOLUTION. IINPUI42 

C I INPU143 

C«««» ONE SET OF K CARDS (8F10.6), K AS ABOVE INPU144 

C i INPUI45 

C XCALB(I), 1= 1,N ... IINPU146 

C X COORDINATE IN CALIBRATION SOLUTION. IINPUI47 

C IINPUI48 

C««H»K ONE SET OF K CARDS (8F10.6), K AS ABOVE m«««n»k««mmi»»«ih«««kk«««m«(hIimPU 149 
C jiNPUISO 

C YCALB(I), 1=1, N ... IINPU151 

C DEPENDENT VARIABLE IN CALIBRATION SOLUTION.. IINPU152 

C IINPUI53 

CNNKW NCASE SETS OF CAROS, EACH SET COMPRISED AS FOLLOWS* INPU154 

C tiNPUISS 

c FIRST CARD (8F10.6I INPU156 

C IINPU157 

C M2 ONCOMING MACH NUMBER IN SOLUTION TO BE IINPUI58 

C COMPUTED. IINPUI59 

C IINPUI 6 O 

C Q2 VALUE OF PERTURBATION PARAMETER IN SOLUTION TOIINPUI 6 I 

C BE COMPUTED. IINPU162 

C IINPU163 

C «•**** ONE SET OF K CARDS (6F10.6), K AS ABOVE INPU164 

C IINPUI65 

C XCHEKdl, 1= 1,N ... IINPUI 66 

C X COORDINATE IN COMPARISON SOLUTION. IINPU167 

C |INPU16a 

C «*•*•« ONE SET OF K CARDS (8F10.6), K AS ABOVE INPU169 

C I1NPU170 



I 


C TCHEK(I>, 1=1 .N ... 

C DEPENDENT VARIABLE IN COMPARISON SOLUTION. 

C 

C «•««»»«» TNE LATTER TWO SETS OF K CARDS ARE 

C • NOTE <• OMITTED WHEN LCHEK = 0 (NO COMPARISON 

C SOLUTION AVAILABLE). 

C 

DIMENSION LOCO(6 1 , LOCI (6 I , LSELCT(6 ) ,TITLE( 20 ) 

DIMENSION XBASE! 200 ) ,XCALB( 200 I , XPERT! 200 ) ,XCKEK( 200 I . 
y. YBASE! 200 ) , YCALB! 200 ) , YPERT! 200 ) , YCHEK! 200 ) 

REAL HC,m,H2 
INTEGER V2 NAME 

COMMON /PARAM/ TITLE, LOCO, LOCI ,L5ELCT,N,NCASE,LSPEC,LECHO,LUNIT, 
y. LCMEK,LPERT,NSELCT,A,B,NAME 


COMMON /PERT/ M0.M1 ,M2,(J0,(}I ,H2,YCR0,tCRI ,YCR2 

COMMON /XY/ XBASE, XCALB. XPERT, XCHEK, YBASE, YCALB, YPERT, YCHEK 

GO TO (100,200,300), ICALL 

100 READ (5,1000) N,NCASE,LSPEC,LECHO,LUNIT, LCHEK, LPERT 

IF (LSPEC .EQ. 0) READ 15,1000) NSELCT,( LSELCT! I ) , 1=1 .NSELCT ) 

IF ILS'EC .EQ. I) READ (5,1000) NSELCT, ( LOCO! I ) ,1=1 .NSELCT 1 , 
y. (LOCH I), 1=1 .NSELCT) 

READ (5,1050) TITLE 

READ (5,1100) NAME 

READ (S.I200) A.B 

READ (5,1200) M0,Q0 

READ (5,1200) (XBASE! I), 1=1 ,N) 

READ (5,1200) ( YBASE! 1 1 ,1 = 1 ,N) 

READ (5,1200) Ml ,Q1 

READ (5,1200) (XCALB( I ),I=1 ,N) 

READ (5,1200) ( YCALB( I ) ,1=1 ,N) 

RETURN 

200 READ (5,1200) H2.Q2 
RETURN 

300 READ (5,1200) (XCHEK( I ) ,1=1 ,N) 

READ (5,1200) ( YCHEKd ),I=1 ,NI 
RETURN 

1000 FORMAT (1615) 

1050 FORMAT (20A4) 

1100 FORMAT IA2) 

1200 FORMAT IBFIO.O) 

END 

SUBROUTINE ECHINP 

DIMENSION LOCO( 6 ) , LOC 1 ( 6 I . LSELCT! 6 ) , TITLE( 20 ) 

DIMENSION XBASE! 200 ) ,XCALB( 200 I , XPERT) 200 ) ,XCHEK( 200 ) , 
y. YBASE! 200 ) , YCALB( 200 ) , YPERT( 200 ) , YCHEK( 200 ) 

REAL MO, Ml, M2 
INTEGERN2 NAME 

COMMON /PARAM/ TITLE, LOCO, LOCI , LSELCT, N.NCASE, LSPEC, LECHO, LUNIT, 
y. LCHEK. LPERT, NSELCT, A.B, NAME 

COMMON /PERT/ MO.MI ,M2.Q0,Q1 ,Q2,YCR0,YCR1 ,YCR2 
COMMON /XY/ XBASE, XCALB, XPERT, XCHEK, YBASE, YCALB, YPERT, YCHEK 
WRITE (6,t<i00) 

WRITE (6,1500) N.NCASE, LSPEC, LECHO, LUNIT, LCHEK, LPERT 
IF (LSPEC .EQ. 0) WRITE (6,1500) NSELCT, ( LSELCT( I ) ,I=t .NSELCT) 

IF (LSPEC .EQ. I) WRITE (6,1500) NSELCT. I LOCO( 1 1 , 1=1 .NSELCT I , 
y. ( LOCI (1). 1 = 1 .NSELCT) 

WRITE (6,1550) TITLE 
WRITE (6,1600) NAME 
WRITE (6,1700) A.B 
WRITE (6,1700) MO,QO 


INPU171 
INPUI72 
INPUI73 
INPUI74 
INPU17S 
INPU176 
INPU177 
IHPU1 7B 
INPU179 
INPU1A0 
INPUIS1 
I))PU162 
INPU1B3 
INPU1B4 
INPU1B5 
INPU1B6 
IMPU1B7 
INPU188 
IHPUiaO 
INPU190 
INPU191 
INPU192 
INPU193 
INPU194 
INPU195 
INPU196 
INPU197 
INPU198 
INPU199 
INPU200 
INPU201 
INPU202 
1NPU203 
INPU209 
INPU205 
INPU206 
INPU207 
INPU20B 
INPU209 
INPU2I0 
INPU211 
ECHIOOI 
ECHI002 
ECHI003 
ECHI009 
ECHI005 
ECHI006 
ECHI007 
ECHIOOa 
ECHI009 
ECHIOlO 
ECHI01 1 
ECHI012 
6CHI013 
ECHI0I4 
ECHI015 
ECHI016 
ECHI017 
ECHIOia 
ECHI019 


<Ti 

LO 


WRITE (6,1700) ( XBASE! I), 1=1 ,N) 

WRITE (6,1700) (YBASE(II,I = 1,N) 

WRITE (6,1700) Ml ,Q1 

WRITE (6,1700) (XCALB(I),I=1,N) 

WRITE (6,1700) (YCALB(I),I=t,N) 

RETURN 

1400 FORMAT ( 1H1 ,25MH«>)/ 

y. IX,1H«,1X,21HLISTING OF INPUT DECK,1X,1H«/ 

y, 1X,25(1H*»)///1 

1500 FORMAT (IX, 1615) 

1550 FORMAT (1X,20A4) 

1600 FORMAT (IX,A2) 

1700 FORMAT (lX,aFl0.6) 

END 

SUBROUTINE banner 
WRITE (6,1300) 

WRITE 16,1310) 

1300 FORMAT ( IHI , I0(/),49X,55(IH* )/49X, 1K«,53X, IH«/ 

y. 49X,1H»,19X,I5HPR0GRAM PERTURB, 19X,1H*/49X,1H*,53X,1H»/ 
y. 49X,tH<*,8X,37HCALCULATES NONLINEAR SINGLE-PARAMETER, 

X 8X,1H*/49X,1H*,53X,1H*/ 

y. 49X,1HA,13X,27HC0NTINU0US OR DISCONTINUOUS, 
y. 13X,1H«/49X,1K>*,53X,tHM/ 

7. 49X,IH«,15X,22HPERTURBATION SOLUTIONS, 16X,1H«/49X,1HM,53X,1HM/ 

■/. 49X,1H»,9X,34HWHICH REPRESENT A CHANGE IN EITHER, 

V. 10X,1HK/49X,1H«,53X,1H«/ 

y. 49X,1H»,13X,27HGE0METRY OR FLOW CONDITIONS, 

V. 13X,1H»/49X,1H«,53X,1H«/ 

y. 49X,1H*,4X,44HBY EMPLOYING A STRAINED -COORDINATE PROCEDURE, 

Y. 5X,1H«/49X,1H»,53X,1H«/ 

y. 49X,1H«,4X,45HUTILIZING A UNIT PERTURBATION DETERMINED FROM, 
y. 4X,1H»/49X,1H*t.53X,1H««/ 

y. 49X,1H»,14X,25HTH0 PREVIOUSLY CALCULATED, 
y. 14X,1H»/49X,1H»,53X,1H«) 

1310 FORMAT I 49X,1H».9X.34H'BASE' AND ‘CALIBRATION’ SOLUTIONS, 

Y. 10X,1H»/49X,1H«,53X,1H»/ 

y. 49X,1H'».4X,45HD1SPLACED FROM ONE ANOTHER BY SOME REASONABLE, 
y. 4X,1H«/49X,1H<*.53X,IKM/ 

V. 49X,1H»|8X,36HCHANGE IN GEOMETRY OR FLOW CONDITION, 
y. 9X,tHft/49X,1H*,53X,1HV 

y. 49X,1H»,53X,1H»/ 

7. 49X.1H'<,2IX,iaHHRITTEN BY.22X, 1H*/49X, 1H»,53X, 1H»/ 

V. 49X,1H«,7X,39H JAMES P. ELLIOTT AND STEPHEN S. STAHARA, 
y. 7X,1H)*/49X,1H»*,53X,1H«/ 

y. 49X,1H*.53X,IH»/ 

y. 49X,1H»,7X,36HNIELSEN ENGINEERING AND RESEARCH, INC., 
y. eX,1HH/49X,1H»,53X,1H*/ 

y. 49X,IH«.I4X,25HM0UNTAIN VIEW, CALIFORNIA, I4X, 1H»/49X, IH».53X, 1H«/ 
y. 49X,55tlH*)) 

RETURN 

END 

SUBROUTINE SCALE (N,X,M,A,B) 

C 

C ENTRY WITH M = 1 CONVERTS FROM PHYSICAL X (0 TO -A ON LOWER 

C SURFAC1 , 0 TO B ON UPPER SURFACE) TO NORMALIZED X (0 < X < 1). 

C ENTRY WITH M=2 REVERSES THE PROCESS. NZ (DETERMINED WHEN M=1 ) 

C CORRESPONDS TO POINT AT NOSE OF BLADE OR AIRFOIL. 


COMMON /FLOREV/ NZ 
DIMENSION X(200) 



•>J 

o 


c 

c. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


GO TO 

1 CONTINUE 
NZ=0 

SO 2 I=2iN 

IF IX(I> .LT. X(I-» n NZ=I 

2 CONTINJE 
00 5 1=1 .N 

IF II .LE. NZI T=-X(I) 

IF tl .GT. NZ) T=Xm 
Xm=(T-A)/IB-A) 

5 CONTINUE 
RETURN 

6 DO 7 1=1. N 
X(I)=ADS((B-A)«X(I)+A) 

7 CONTINUE 
RETURN 
END 

SUBROUTINE LOCATE (N,X> Y. YCRIT.IGRAD, LNIN, LMAX.NCRIT. LCRIT.XLOC ) 


.OPERATES ON THE INPUT ARRAY Y, LOCATING NINIHUM AND MAXIMUH I 
VALUES, AND ALL CRITICAL POINTS (Y=YCR1T) FOR WHICH DY/OX 1 IN I 
PHYSICAL COORDINATES) HAS ALGEBRAIC SIGN GIVEN BY IGRAD. NCRIT I 
IS NUMBER OF CRITICAL POINTS. POINTS FOUND ARE STORED IN THE ARRAYl 
XLOC AS FOLLOWS! I 


XLOCI 1 ) = MINIMUM PT. 
XLOC(Z) * MAXIMUM PT. 
XL0C13) = CRITICAL PT, tl 

XLOCI 6) = CRITICAL PT. t4 


DIMENSION XI 200 ) . Y( 200 ) . LCRITI 4 ) .XCRITI ) .XLOCI 6 ) 

COMMON /FLOHEV/ IREV 

IFL0H=-1 

LMIN=1 

LMAX=1 

ISTART=2 

IF (IREV .EG. 0) GO TO 5 
LMIN=2 
LHAX=2 
ISTART=3 
5 CONTINUE 
NCRIT=0 

DO 100 I=ISTART.N 

IF (IREV .NE. 0 .AND. I .EQ. N) 'GO TO 10 
IF (YIII .GT. Y(LMAXI) LMAX=I 
IF (Yin .LT. YILMIN)) LMIN=I 
10 CONTINUE 

IF (Itm .GT. YCRIT .ANO. Y(I-1) .GT. YCRIT) .OR. 

X IY(I) .LT. YCRIT .AND. YII-D .LT. YCRIT)) GO TO 100 
IF (I .GT. IREV) IFLOW=l 

IF UY(I)~Y(I~1 ))KFLOAT(IFLOU»IGRAD) .LT. 0.0) GO TO 100 

NCRIT=NCRIT*t 

LCRIT(NCRIT)=I-1 

SLOPE=(X(I)-X(I-t ) )/(Y(I)-Y(I-l )) 

XCRITI NCRIT )=X( 1-1 )+SLOPE»I YCRIT-YI I-l) ) 

100 CONTINUE 

XLOCd )=X( LMIN) 

XLOC(2)=X(LMAX) 

IF (NCRIT .EG. 0) RETURN 


C 

C 

C 

C 


c 

c 

c 

c 


c 

c 

c 


DO 200 1 = 1, NCRIT I LOCAOA^i 

200 XLOC(I*?)=XCRIT(I) I LOCAO^iS 

RETURN I LOCAO«i6 

end I LOCAO!*7 

SUBROUTINE SORT (N.X.ISEQ) ISORTOOl 

ISORT002 

....ARRANGES THE SET X( 1 ) , X( 2 ) , ... , X(N) IN A MONOTONE INCREASING 1S0RT003 
SEQUENCE. ISEQ GIVES ORDER OF SUBSCRIPTS IN REARRANGED SEQUENCE. I S0RT00<( 

ISORTOOS 

DIMENSION X(B),ISEQ(8) ISORT006 

NM1=N-1 IS0RT007 

DO 1 1=1 ,N I SORTOOa 

1 ISEQ(I) = I ISORT009 

10 ITEST=0 iSORTOlO 

DO 100 1=1. NMt ISORTOII 

IF (X(I) .LE. X(I+1)) GO TO 100 ISORT012 

XSAVE=X(I) ISORTOII 

X( I )=X 1+1 ) ISORT019 

X(I + 1)=XSAVE ISORTOtS 

ISAVE=ISEQ(I) IS0RT016 

ISEQ( I l=ISEQ(I+l ) ISORT017 

ISEQ( 1 + 1 )=ISAVE l.SORTOia 

ITEST=1 ISORT019 

100 CONTINUE ISORT020 

IF (ITEST .EQ. 1) GO TO 10 1S0RT021 

RETURN ISORT022 

end ISORT023 

SUBROUTINE INTERP (N.X.Y.XI.YI ) lINTEOOl 

IINTE002 

GIVEN THE SET OF POINTS X(I), Y(I), 1 = 1 ,N, ANO THE SET XKJ), 1 1NTE003 

J=1,N, USES LINEAR INTERPOLATION TO COMPUTE THE SET VKJ), J=1 ,M. IINTE00<* 

IINTE005 

DIMENSION X(200),Y(200).XI(200),YI(200I IINTE006 

HM1=N-1 I1NTE007 

JSTART=1 lINTEOOa 

DO 100 1=1 ,N I 1NTE009 

IF (XIII) .LE. X(1l) GO TO 10 IlNTEOlO 

IF (XIII) .GE. X(N)) GO TO 20 lINTEOIl 

GO TO 30 I1NTE012 

10 J=1 IINTE013 

GO TO 95 IINTE0I9 

20 J=N-1 IINTE015 

GO TO 95 IINTE016 

30 CONTINUE IINTE017 

DO 90 J=JSTART,NMI IINTE018 

IF (XKI) .NE. XUn GO TO 90 IINTE019 

YIII) = YIJ) 1INTE020 

GO TO 1 00 I INTE02I 

90 IF (XKI) .GT. X(J) .AND. XKI) .LT. XIJ + D) GO TO 95 IINTE022 

90 CONTINUE I INTEC23 

95 SLOPE = (Y(J + 1 )-Y(J))/(X(J + ) )-X(J)) I INTE029 

YI( I )=Y( J )+SLOPE!(( XI( I )-X( J ) ) I INTE025 

JSTART=J IINTE026 

100 CONTINUE IINTE027 

RETURN IINTE028 

END IINTE029 

SUBROUTINE STRAIN I N ,NSEG,XFIX ,XIN . PARM.XOUT ) ISTRAOOl 

ISTRA002 


COMPUTES STRAINED COORDINATE FROM INPUT ARRAY XIN, USING PIECEWISE I STRA003 
LINEAR STRAINING WITH NSEG LINEAR SEGMENTS. FOR UNIT STRAINING, ISTRA009 



C INPUT VALUE OF PARM IS 1.0; FOR GENERAL CASE, 

C PARtI = tt)2-Q0)/(qi-Q01. 

C 

DIMENSION XFIX(6l,XIN(200),XOUTt200 ) 

COMMON /COEFF/ C(7),D(7) 

JSTART=1 

00 50 1=1 .N 

00 ')0 J=JSTART,NSEG 

IF IXINdl .GE. XFIX1J) .AND. XIN( 1 1 ,LE. XFIXU+ll) GO TO 45 
40 CONTINUE 

45 XOUTri)=XINm + PARM«(C(JJ + tDIJ)-t.01»XINIin 
JSTART=J 
SO CONTINUE 
RETURN 
END 

SUBROUTINE MONO (N,L,X.Y) 

C 

C CHECKS POINTS IN VICINITY OF A CRITICAL POINT FDR MONOTONE 

C BEHAVIOR, AND ADJUSTS VALUES TO GIVE A LINEAR PROFILE. 

C 

DIMENSION L(4).X(ZOO),Y(20O) 

DO ICO 1=1, N 
LS=L( 1 1 
Y1=Y( LS-1 ) 

Y2=Y1 LS) 

YJ=Y( LSH ) 

Y4=Y( LS»2 ) 

IF tlYi .LT. Y2) .AND. ( Y2 ,LT. YJ) .AND. IY3 .LT. Y4 1) GO TO 100 

IF KYI .GT. Y2) .AND. 1 Y2 .6T. Y3) .AND. IY3 .GT. Y4 1) GO TO 100 

X1=X( LS-1 » 

X2=X( LSI 
X3=X(LSe1 » 

X4=X( LS+21 

SL0PE=(Y4-Y1 )/(X4-X1 ) 

Y( LS)=YKSL0P£»(X2-X1 > 

YtLS+1 1=Y1+SLOPE*IX3-X1 ) 

100 CONTINUE 
RETURN 
END 

SUBROUTINE PLOT (N.LPERT) 

C 

C CREATES FILE OF COMMANDS FOR PROGRAM 'TQPDRAW AT STANFORD CENTER 

C FOR INFORMATION PROCESSING (S.C.I.P.1. CALLED ONLY ONCE IN MAIN 

C PROGRAM AND MAY BE DELETED OR REPLACED. 

C 

DIMENSION XO(200),X1 (200 ),X2( 200),X3( 200), 

X YO( 200),Y1(200J,Y2(200),Y3I200) 

COMMON ./PERT/ MO, Ml ,M2,Q0,O1 , 02 , TCRO , YCR1 ,YCR2 
COMMON 'Xy/ X0,X1 ,X2,X3,Y0,Y1 ,Y2,Y3 
LOGICALOl SYM(3) /I HT , 1 HA , 1 HM/, SUBI3) /1HG,1HG,1H /, 
y. SO /1HB/. SI /1HC/ 

REAL MJ, Ml, M2 

DATA lUALL /O/, XL,XR0,XT0,XR1 ,XT1 ,XR2,XT2 
7. /O.O, 0.12, 0.14, 0.16, 0.1B, 0.2, 0.24/ 

IF (ICALL .EQ. 0] WRITE 14,10001 
ICALL=1 

CALL LIMITS IN,YMIN,YMAX) 

. WRITE (4,1050) 

WRITE (4,11001 YMIN,YHAX 
WRITE 14,12001 


SIRAO05 
S1RA006 
5TRA007 
STRAOOa 
STRA009 
SIRAOlO 
SIRA01 1 
S1RA01 2 
SIRA01 3 
STRA0I4 
STRA015 
S1RA016 
STRA01 7 
STRAOia 
STRA019 
MQN0001 
M0MO002 
MDN0003 
MQN0004 
MOMO005 
MOMO006 
MOM0007 
MOM0009 
MOMO009 
MOM0010 
MOMO01 I 
HOMO012 
MONO015 
MOMO014 
MOMO015 
HON0016 
MOM0017 
MOMOOia 
MON0019 
HONDO 20 
MOMO021 
MOMO022 
MONO023 
MONO024 
PLOTOOl 
PLOT002 
PLOT003 
PLOT004 
PLOT005 
PLOT006 
PLOT007 
PL01008 
PLOT009 
PLOT010 
PLQT01 1 
PL0T012 
PLOT013 
PLOT014 
PLOT01S 
PLOT016 
PLOT01 7 
PLOT018 
PLOT019 
PLOT020 
PLOT021 


WRITE 

(4, 

1 100 1 

WR] TE 

(4, 

1 350) 

WRI (E 

(4, 

(3b0J 

WRITE 

(4, 

1 370 ) 

WRITE 

(4, 

1400) 

IF ( LPERT 

.NE . 

IF ( LPERT 

■ Eq. 

IF ( LPERT 

.ME. 

WRITE 

(4, 

5000 1 

WRITE 

(4, 

1 700 ) 

WRITE 

(4 , 

1530 1 

WRITE 

(4, 

5000 ) 

WRITE 

(4, 

1600 1 

WRITE 

(4, 

1530) 


EO. 31 WRITE (4,1505) 02 


to CONTINUE 

WRITE (4,50001 XL , YCR2 , XR2 , YCR2 

WRITE (4,15101 

WRITE (4,15501 XT2,YCR2 

WRITE 14.50001 1 X01 1 ) , YOl I ) , 1 = 1 ,N ) 

WRITE 14,1 7001 

WRITE (4,50001 ( XI ( I ) ,Y1 ( 1 1 , 1=1 ,N I 
WRITE 14,16001 

WRITE 14,50001 ( X2( I ) , Y2( I ) ,1=1 ,N ) 

WRITE (4,19001 

WRITE (4,5000 ) ( X31 1 ) , Y3( I ) , 1 = 1 ,N ) 

WRITE (4,15(01 
RETURN 

1000 FORMAT 1 45H//PERTPLOT JOB JIM ELLIOTT ' ,REGI0N=51 2K/ 

•/. 19H//0RAW EXEC TOPDRAW/ 

X 40HSET DEVICE VERSATEC CONTINUOUS INTENSITY/ 

■/. 18HSET CARD LENGTH 80/ 

Y. 15HSET FONT DUPLEX) 

(050 FORMAT 1 9HNEW FRAME/ 

X 27H5ET TICKS TOP OFF RIGHT OFF/ 

X 22H5ET WINDOW X 2 7 Y 2 6/ 

X 2011SET SYMBOL 6P SIZE 1 ) 

1100 FORMAT (10HSET LIMITS X 0 1 Y,2F5.1) 

1200 FORMAT (51HTITLE 4.5 9.5 CENTER SIZE 3 SPACES 7 'PLOT OF C2P3V 
X 19HCASE ' LLL LL CLC ) 

1300 FORMAT (39HTITLE 4.5 8.7 CENTER SIZE 1.5 SPACES 38, IX, 

X 24H'FULL F2HOLST3 NONL. ') 

1350 FORMAT 1 4HMORE , IX, 1H ' , A1 , 1X,5HPERT. , 3X,A1 ,5H2B3 =,F6. 3,3X,A1 , 

X 5H2C3 =,F6.3,1H' ) 

1360 FORMAT ( 4HCASE , 1 X, 24H ' LLL GC LLLLC LLL ') 

1370 FORMAT ( 4HCASE , 1 X , 1 H ' ,At ,2X, 3HLLL,4X,A1 ,3HCLC, I tX,A( ,3HCLC,8X, 1H’ ) 
1400 FORMAT (25MTITLE 0.8 5 SIZE 2 'C2P3'/ 

X 4HCASE,1X,6H' CLC’/ 

X 23HTITLE BOTTOM SIZE 2 'X' 1 

1500 FORMAT 147HTITLE 4.5 0.5 CENTER SIZE 1.5 SPACES 17 'M203 =, 

X F6.3,3X,A1 ,2H =,F6.3,1H'/ 

X lOdCASE ' CSC,l IX, 1HG,6X,1H' ) 

1505 FORMAT (34HTITLE 4.5 0.5 CENTER SIZE 1.5 'M =,F6.3,1H') 

1510 FORMAT (6HJ0IN 1) 

1530 FORMAT ( 5HTITLE , 1 X , 2( F6 . 3 , 1 X ) , 

X 27HDATA SIZE 1.0 ' 1 C42P350X1 )2 , A1 , 2H3 ' / 

X 4HCASE,1X,16H' CCLCCC C CLC'l 

(550 FORMAT ( 5HTITLE , IX , 2( F6 . 3, 1 X ) ,25HDATA SIZE 1.5 'C4ZP350»*r/ 

X 16HCASE ' CCLCCC C ) 

1700 FORMAT 115H5ET INTENSITY 3/ 

X 11HJOIN 1 DOTS/ 
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y. 15H5ET INTENSITY 2) 

1800 FORMAT M3HJOIN 1 DASHES) 

1900 FORMAT 19HPL0T) 

5000 FORMAT ( 51 2( F6 . 3 . 1 X) , 1 H ; ) ) 

END 

SUBROUTINE LIMITS ( N, YMIN, YMAX ) 

C 

C SEARCHES FOUR DATA ARRAYS Y0,Y1,Y2,Y3 FOR MINIMUM AND MAXIMUM. 

C CALLED ONLY BY PLOT SUBROUTINE. 

C 

DIMENSION X0( 200 I ,X1 ( 200 ) ,X2( 200 ) ,X3< 200 ) . 

•/. Y0( 200 ) ,YU 200 ) ,Y2( 200 ) >Y3< 200 ) 

DIMENSION Z(800) 

COMMON /XY/ X0.X1 ,X2,X3,Y0,Y1 .Y2.Y3 
EQUIVALENCE ( Y0( 1).Z( 1 I) 

YMIN=Z' I ) 

YMAX=Z 1 ) 

DO 10 1=1 .A 

JSTART = 200»«( 1-1 ) + 1 

JSTOP=JSTARTtN-1 

DO 10 J=JSTART,JSTOP 

IF (Z(J).GT.YMAX) YMAX=Z(J) 

IF (Z(J).LT.YMIN) YMIN=Z(J) 

10 CONTINUE 
YSAVE=YMAX 
YMAX=YMIN 
YMIN=YSAVE 
CALL ROUND ( YMIN ) 

CALL ROUND ( YMAX ) 

RETURN 

END 

SUBROUTINE ROUND (Y) 

C 

C ROUNDS Y LIMITS FOR OUTPUT IN F5.1 FORMAT. CALLED ONLY BY 

C SUBROUTINE LIMITS. 

C 

Z=ABS( Y) 

IF ( 10.»Z-INT( 10.»Z).LT. .5) Z=Z+ . 05 
IF (Y.GT.O. ) GO TO 1 
Y=-Z 
RETURN 
1 Y=Z 
RETURN 
END 

SUBROUTINE UPLOW ( A,B,XIN,K,N,XOUT .FLAG ) 

C 

C CONVERTS NORMALIZED ARRAY XIN TO PHYSICAL ARRAY XOUT AND FLAGS 

C POINTS OH LOWER SURFACE WITH A •«’. 

C 

DIMENSION XIN(K),X0UT(8) 

LOGICALHI FLAG(8), BLANK/1 H /. STAR/I H»/ 

XN05E=-A/(B-A) 

DO 1 1=1 iN 
FLA6( I )=BLAHK 

IF (XIN(I) .LT. XN05E) FLAG(I)=5TAR 
XOUT( I ) = AB5( (B-A)»XIN( I ) + A) 

1 CONTINUE 
RETURN 
END 
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APPENDIX C 


C 

H 

i 


n 


N 


q 

q 

q 

Q 

Q 

Q 

Q 


c 

o 


c 


o 


p 


(s , t) 


LIST OF SYMBOLS 


blade chord, m 

blade spacing for nonstaggered cascades, m 

invariant point index; eq. (6) ; also, index for 
Lagrangian coefficients; eq. (22) 

dummy index; eq. (20) 

two-dimensional full potential operator; eq. (2) 

linear operator representing first-order perturbation 
of two-dimensional full potential equation; eq. (4) 

linear operator representing first-order perturbation 
terms arising from coordinate straining; eq, (9) 

Lagrangian coefficients; eq. (22) 

total number of shock points and high-gradient 
maxima points; eq. (24) 

total number of invariant points, equal to n+2; 
eq. (24) 

arbitrary geometric or flow parameter to be perturbed; 
eq. (13) 

calibration flow value of q; eq. (9) 
base flow value of q; eq. (3) 

approximate flow solution for arbitrary flow quantity; 
eq. (1) 

calibration flow solution for value q of arbitrary 
parameter; eq. (8) 

base flow solution for value q of arbitrary parameter; 
eq. (1) ^ 

linearized perturbation solution per unit change of 
perturbed parameter; eq. (1) 

strained (x,y) coordinates; eq. (5) 
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(x,y) 


nondimensional blade-“fixed orthogonal coordinates; 
eq. (5) , normalized by. C 

(x,y) nondimensional blade-fixed orthogonal coordinates 

related to calibration solution; eq. (9) 

(x, ,y ) straining functions associated with (x,y) coordinates; 

eq. (6) 

(Xf ,y^ ) straining functions associated with ith invariant 
■^i "^i point; eq. (6) 

(6x.,6y.) unit displacements in (x,y) directions associated with 
ith invariant point; eq. (6) 

c 

6x. unit displacement in x direction between base and 

calibration flows of the ith invariant point; eq. (18) 

€ perturbation change of geometric or flow parameter; 

eq. (17) 

perturbation of geometric or flow parameter between 
base and calibration flows; eq. (18) 

^ nondimensional total velocity potential; eq. (2), 

normalized by CV^ 

$ nondimensional base flow velocity ptential; eq, (3), 

normalized by CV 

-I CO 

nondimensional perturbation velocity potential; eq, (3), 
normalized by CV 

00 

Subscripts 

0 denotes base flow quantities 

1 denotes perturbation quantities 

c denotes quantities associated with calibration 

flow 
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Figure 2.- Summary of various flows and straining 
functions considered 
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Figure 3.- Comparison of perturbation (0) and non 
linear { — ) surface pressures for a thickness- 
ratio perturbation of a nonlifting cascade 
of biconvex profiles with H/C = 1.0 
at M =0.80 
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Figure 5.- Comparison of perturbation (O) and non 
linear ( — ) surface pressures for a thickness- 
ratio perturbation for an isolated NACA OOXX 
airfoil at = 0.820 and a = 0° for 
extreme solution extrapolation 
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Figure 7.- Comparison of perturbation (O) and 
nonlinear ( — ) surface pressures for an 
angle-of-attack perturbation of an 
isolated NACA 0012 airfoil at 
= 0.70 
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Figure 10.- Comparison of nonlinear ( — ) sur- 
face pressures with perturbation results 
using quadratic (O) and linear piece- 
wise-continuous (*) straining functions 
for a thickness-ratio perturbation of 
an isolated NACA OOXX airfoil at 

M = 0.820 and a = 0° 
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Figure 11. ~ Comparison of nonlinear ( — ) sur- 
face pressures with perturbation results 
using cubic (O) and linear piecewise- 
continuous (*) straining functions for 
an angle-of-attack perturbation of an 
isolated NACA 0012 airfoil at 
M =0.70 
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Figure 12.- Comparison of nonlinear { — ) sur- 
face pressures with perturbation results 
using quadratic (O) and linear piece- 
wise-continuous (*) straining functions 
for a Mach number perturbation of an 
isolated NACA 0012 airfoil at 
a = 0° 
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Figure 13.-* Comparison of nonlinear ( — ) sur*- 
face pressures with perturbation results 
using cubic (O) and linear piecewise- 
continuous (*) straining functions for 
an angle-of-attack perturbation of an 
isolated NACA 0012 airfoil at 

M =0.70 
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